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varies over its range 
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directed line segment between two connected points on the 
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reference trajectory position, velocity, acceleration time- 
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attachment points and spreader-bar endpoints) (see fig. 1.1) 

transformation of vector coordinates from reference frame "b" 
to frame "a"; the subscripts can be any of H, N, t for the 
frames used here 
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helicopter thrust limits 
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max 
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components of in level-heading axes 
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weight; subscript (t, b, 1,2) indicates body 
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®b»®h 


aerodynamic angles defining direction of air velocity vector in 
body axes: 

Va = Va(cos a cos B i + sin 8 j. + sin o sin 8 k) ; subscripts 
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(e,x) direction angles of a vector in triangle axes: 
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tion with F34 = -A 




Pb>PT*pTx 
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direction angles (angle from vertical and heading) of ut in 
level heading axes: 
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angle defined in context (sec. 6) 
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magnitude 
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Euler angles of transformation from level-north axes to body 
axes; subscripts indicate the body 
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SUMMARY 


The characteristics of a dual-lift helicopter system in equilibrium flight 
along any reference trajectory are examined. The system consists of the cargo 
attached by cables to the endpoints of a spreader bar which is suspended by tether 
cables from the two helicopters. Results are given for the orientation angles of 
the suspension and its internal forces, and for the helicopter thrust vector 
requirements under general circumstances, including nonidentical helicopters, and 
any static or accelerating equilibrium flight condition. These results provide a 
basis for coordinating the system in flight. 

Analytical results are derived from the force balance equations at the centers 
of gravity of the load and helicopter, and at the endpoints of the spreader bar. 

The analysis does not consider the attitude degrees of freedom of the load and 
helicopters, but assumes that these are stable and that the aerodynamic forces in 
equilibrium flight on these bodies can be calculated Independently and to sufficient 
accuracy from aerodynamic models or that they can be measured in flight. This 
assumption permits closed-form solution of the equations with the equilibrium aero- 
dynamic forces appearing parametrically. The nature and magnitude of these forces 
are examined for example helicopters and loads. 

It is found that there are arbitrarily many equilibrium orientations of the 
system for any one situation (any given helicopter pair, load, and flight condition) 
with three angles which can be selected by the pilot or autopilot. These can be 
taken as (1) the spreader-bar heading relative to the ground-track direction (forma- 
tion angle), (2) the spreader-bar tilt relative to the apparent gravity of the sus- 
pended load modified to include the load aerodynamic specific force (force per unit 
mass), and (3) the tilt of one of the tethers relative to the spreader bar. It is 
shown that thrust requirements vary strongly with the tether angle but are nearly 
invariant with the other two independent variables; hence the tether angle can be 
selected to minimize the sum of the required thrust magnitudes. This choice maxi- 
mizes the system thrust margin and also minimizes fuel rate in the case of identical 
helicopters. The spreader-bar tilt controls the distribution to the two helicopters 
of the thrust requirements imposed by the suspended cargo. If the helicopters are 
identical then zero tilt results in equal distribution and, more generally, the tilt 
can be selected for any distribution for which no cable collapses and for which both 
helicopters have positive thrust margins. It is found that the thrust requirements 
can always be distributed in the same ratio as the helicopter thrust limits. Last, 
the formation angle can be selected as in conventional formation flying at some 
fixed value away from the longitudinal and transverse arrangements (e.g., at 45°). 
However, with fully automatic flight control, it may be possible to select this 
angle from a larger range or vary it in flight. The present study does not assume a 
uniquely suitable choice of formation angle, but examines its considerable influence 
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on the system's inertial orientation and the variation of this orientation with 
maneuvering along a reference trajectory. 


1 . INTRODUCTION 


The airborne transport of loads slung beneath two or more helicopters has been 
considered periodically since the 1950s, following the success of single helicopter 
slung load operations. This technique would multiply the maximum load-carrying 
capacity available from existing helicopters used singly. Interest in dual-lift is 
represented by a few civil applications of the technique and by studies sponsored by 
the Department of Defense. The use of two helicopters has been a feasible and 
economic alternative to acquiring existing helicopters of larger payload in some 
commercial operations with simple load suspension requirements (described in 
ref. 1), and has been advocated by Carter et al. (refs. 1 and 2) as an alternative 
to developing new helicopters with payloads exceeding those of current helicopters. 

In a 1957 feasibility study of suspension systems to couple two to six helicop- 
ters with a load, VERTOL (ref. 3) recommended a system of spreader bars suspended 
below the helicopters over (1) a system with rigid connection of the helicopters, or 
(2) a system without separation restraint in which the load was attached solely by 
cables to each helicopter. This study was based on tandem rotor helicopters each 
with a maximum payload of 2 tons. A decade later, Sikorsky Aircraft studied the 
feasibility of twin helicopter systems (refs. 1 and 4) and also recommended the 
spreader-bar suspension over several other methods of suspension. This system 
(fig. 1.1(a)) consists of a spreader bar hung below the helicopters by tether cables 
with the load attached by cables to the endpoints of the spreader-bar. The dimen- 
sions of the system are chosen to minimize spreader-bar weight while ensuring a safe 
helicopter separation; this results in proposed spreader bar lengths of 1.25 to 
2 rotor diameters, depending on the control system, and in bridle cable angles of 
45° to 60°. This study was based on the single-rotor CH-54B with a 10-ton payload, 
and the recommended system was subsequently demonstrated in flight tests which were 
confined to low speeds and within the test airfield (refs. 1 and 5). 

There are also several examples of civil twin-lift operations; these have been 
flown over longer distances and speeds to 60 knots, but have generally been limited 
to loads of sufficient length (poles, logs, pipes, cables, etc.) to be suspended 
just below the helicopters without using a spreader bar (fig. 1.1(b)). This 
arrangement is inherently more stable than the system of figure 1.1(a) which accom- 
modates more typical loads that are significantly shorter than the minimum spreader 
bar length. At present, interest in further development of dual-lift operations 
appears in proposed extensions of the flight tests reported in reference 5. Thus, 
dual-lift has been considered at several stages in the expansion of single helicop- 
ter payloads, but has yet to be developed into a proven operational system. A key 
step in achieving this goal is development and demonstration of a suitable control 
system, which is made feasible, in part, by today's digital flight control 
technology . 
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HELICOPTER #2 




(b) LONG LOADS 

Figure 1.1.- Dual-lift configuration. 


The flight demonstration (ref. 5) tested the system at hover and at speeds to 
20 knots and pilots encountered high work load, particularly during maneuvering when 
the spreader-bar tilt must be varied as necessary to maintain system equilibrium 
with invariant load distribution to the two helicopters. This work load indicated a 
potential need for coordination aids and for automatic stabilization and control to 
avoid excessive or even unstable transient motion during maneuvering. Indeed, it is 
anticipated in reference 1 that automatic control is necessary if the full opera- 
tional potential of the system is to be achieved, and a master-slave stabilization 
and control system is proposed (refs. 1 and 6). In addition, linear system stabil- 
ity analyses and simulation studies are reported in references 3, and 7 based on 

simplified system and motion models and for hover. 
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The present study was undertaken to explore the characteristics of the system 
shown in figure 1.1(a) in equilibrium flight along any reference trajectory as a 
basis for coordinating the system during maneuvering throughout its flight enve- 
lope. Some basic rules for appropriate tether directions and spreader-bar attitude 
during static equilibrium and turning flight in the case of identical helicopters 
and the effect of bar tilt on load distribution are already recognized in the flight 
tests cited above. The present paper seeks to extend these approximate rules to 
general maneuvering flight and to include the effects of aerodynamic forces on the 
system's rigid bodies. The basic assumptions are discussed in section 2 , the system 
equations and method of solution in section 3> and the aerodynamic forces and c.g. 
accelerations of the system's bodies in equilibrium flight in section 4. Results 
are then given for the bridle cable forces (sec. 5), the inertial attitude of the 
spreader-bar-load triangle (sec. 6), and spreader bar and tether forces and thrust 
vector requirements (sec. 7), along with an analysis of the tether angle optimiza- 
tion (sec. 8), the distribution of the suspended load to the two helicopters 
(sec. 9) and the variation in system equilibrium geometry and forces with tether 
angle (sec. 10). 

Finally, we remark that this study assumes load stability. However, the speed 
regime over which this assumption is satisfied for various loads with the suspension 
system of figure 1.1(a) is unknown; it is likely to be well below the power-limited 
maximum speed of the helicopters in view of the system's lack of load yaw restraint 
and long suspension distance. The inherent load stability at higher speeds was not 
considered in the present design; since this has been the most influential factor in 
the geometry of advanced single helicopter load suspensions, it may also strongly 
influence future development of dual-lift suspensions. Thus, future systems may be 
significantly altered from that of figure 1.1(a) and the system of the present study 
may be limited to low speeds and favorable loads. 
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2. NOTATION AND ASSUMPTIONS 


The dual lift system is assumed to be arranged as shown in figure 1.1(a). It 
consists of four rigid bodies (two helicopters, a spreader bar, and a suspended 
load) linked by tethers and bridle cables. The bridle cables and spreader bar form 
a rigid isosceles triangle with angle 6 in the range of 45° to 60°. The most com- 
monly considered system consists of identical helicopters with equal loading, but 
the analysis allows for unlike helicopters and unequal loading. A truncated system 
using two helicopters and tethers to transport sufficiently long loads (pipes, logs, 
bridge sections, etc.) as shown in figure 1.1(b) is also potentially useful. The 
arrangement in figure 1.1(a) was recommended in early studies (refs. 3 and 4), from 
among several choices, based on the gross properties of feasibility and thrust 
penalties. However, the development of suspensions for single helicopters for 
military cargo transport has been greatly influenced by load restraint and stabil- 
ization considerations, with the object of stabilizing all loads in a minimum-drag 
attitude over the power-limited speed range of the helicopter (refs. 8-10 and their 
reference lists). These objectives are not yet integrated in the dual-lift suspen- 
sion of figure 1.1(a); as a result, future designs may differ from the arrangement 
assumed here. 


An estimate of the suspension system dimensions can be given. The flight 
demonstration (ref. 5) with twin single-rotor helicopters used a system with a width 
(spreader-bar length) twice the rotor diameter (2D) and load-to-helicopter height 
2D; this is expected to be a maximum size. This size can be decreased through the 
use of reliable automatic station keeping to preclude rotor contact; a minimum size 
of 1.25D by 1.25D is predicted in reference 1. Rotor diameters of existing helicop- 
ters fall in the range of 30-125 ft (ref. 11). 


The sizes of the component masses relative to the load mass can be estimated. 
For given helicopters and mission profile let PL1 and PL2 be the payloads of heli- 
copters 1 and 2, respectively, and assume that dual-lift is both necessary (cargo 
weight exceeds 50j of the combined payload) and possible (suspended weight is less 
than the combined payload): 


W, > max(PL1, PL2) > (PL1 + PL2)/2 

Xr 

+ Wb < (PL1 + PL2) 

where Wj^ and Wj^ are the weights of the load and spread-bar, respectively. 
Express W|, as a fraction of the combined payload: 


( 2 . 1 ) 


Wjj = 6jj(PL1 + PL2) (2.2) 

For efficient use of dual-lift, the combined payload is much larger than the bar 
weight; e.g., < 0.1 in reference 1. Combining equations (2.1) and (2.2) gives 
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(2.3) 


where the upper bound corresponds to the lightest load. Therefore, 

if: Bi, << 1 then: m. << m, (2.^0 

b b t 

Thus, the bar weight is an order of magnitude smaller than the cargo weight, pro- 
vided that equation (2.1) and the condition of equation (2.4) are satisfied. Next, 
express the helicopter weights as a fraction of the payload: 


Wi + W2 = Bj^(PL1 + PL2) 


Then, using equation (2.1), 


1 - 8v 


hlh 

W. 


< 2Bu 


(2.5) 


( 2 . 6 ) 


Noting that helicopter gross weights are about twice their empty weight (ref. 11) it 
can be estimated that is in the range of 1 to 3 for external load missions. 
Thus, the sum of helicopter masses is roughly the same as the load mass, provided 
that the previous assumptions hold and that the vehicle payload is above, say, one 
third of the vehicle weight. 


ui ^ -r 

if: B. << 1 and 1 < Bv, < 3 then: 1 < — < 6 (2.7) 

b h m. 

l 

and, for twin vehicles, each mass is in the range of 0.5 mj^ to 3 mj^. The results in 
equations (2.4) and (2.7) indicate the range of relative masses of the component 
bodies of interest for the dual-lift system. 

The dual-lift configuration is defined by the positions of five points on the 

suspension system, which are numbered for this work as shown in figure 1.1(a). The 

notation for cable and spreader-bar directions and forces is based on this enumera- 
tion; i.e., Ri, Ri.1 , Lij, ui.1 , and FiJ refer to the position of the ith point, 
the link between points i and J, and its length and direction vector, and the scalar 
force carried by the link. The (scalar) link force Fij is taken to be positive 
(negative) if the link is in tension (compression); it is applied to the point Ri 
(R.l) in the direction ui.1 ( - ui.1 ) . 

The component masses (suspended load, spreader-bar, and the two helicopters) 
and their accelerations and applied forces are indicated by the symbols t, b, 1, 

and 2 appended to the notation; e.g., m 2 , a2, fa2, and refer to the mass, 

acceleration, aerodynamic specific force, and specific apparent load vector of 
helicopter No. 2. 
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Forces appear in the analysis principally as "specific” force; that is, force 
per unit weight. In addition, for any component body it is useful to define the 
specific "apparent load" f; that is, the sum of the apparent gravity, g - a, and the 
applied external specific forces exclusive of control or thrust forces: 

f = g + fa - a (2.8) 

where fa is the aerodynamic specific force on the body. This is dominated by 
gravity in the present context. For the load, ft is the specific apparent sus- 
pended load supported by the system, and, in general, mf is the apparent load owing 
to the weight, aerodynamic force, and acceleration of a component mass which must be 
supported by the helicopter engines. Indeed, the required thrust vector sum is 
given by the sum of apparent loads: 

II + T2 = -(m^^ft + mj^fb + m^f1_ + m^^fg) 

and most results are given in terms of the apparent loads of the system's component 
bodies. 

Several orthogonal right-handed reference frames are used in this work. Their 
axes are denoted (i, J[, k) with a subscript (N, H, t) to indicate the frame (level- 
north or level-heading axes, or body axes for the triangle). Level-north and level- 
heading axes are both local vertical frames with ijj and i^^ directed north along 
the reference trajectory ground track, respectively. The triangle formed by the 
spreader bar and bridle cables is a rigid body within the assumptions of this 
work. The triangle axes are body axes attached to the triangle (figure 2.1) and 
oriented with i^. along the spreader bar (along R43 ) and k^ perpendicular to the 
spreader bar in the plane of the triangle and pointed toward the load attachment 
point. Most of the analysis is carried out in this reference frame. 

Subscripts are also used to indicate a vector given by its components in a 
reference frame. Thus, a vector v given in reference frame "a" is denoted by 
v^ where 

Va = (v-ia* V.4, v-k^) 

and the subscript "a" can be any of those noted in the preceding paragraph. 

Transformations of vectors between orthogonal right-handed frames "a" and "b" 
are denoted T^ 1 ^; that is, if v is any vector and v^ and V]^ refer to its compo- 
nents in reference frames "a" and "b", then 


Transformation matrices T^ can be given as a sequence of single axis rotations, 
Ej(o), where i indicates ihe axis of rotation (see fig. 2.2). Euler-angle trans- 
formations suffice for the present work and are given by the following rotation 
sequence and customary aeronautics notation for the angles: 
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UNIT VECTORS: 


u = -sin e + cos e sin X + cos e cos X 
= (-sin e, cos e sin X, cos e cos X) 

ANGLE NOTATION FOR UNIT VECTORS OF INTEREST: 


VECTOR 

UNIT VECTOR 

X 

6 

fC 


Xg 


R13 

u13 

^13 

^13 

R24 

u24 

X 24 

CM 

-11 

Hi 

Xi 


-T2 

u 2 

X 2 



Figure 2.1- Triangle body axes and direction angles e,X 

to triangle axes. 


relative 



SINGLE-AXIS 

ROTATIONS 



1 

0 

0 

Ei(a) = 

0 

cos 0 

sin 0 


0 

-sin 0 

cos 0 


cos 0 

0 

-sin 0 

E2(a) = 

0 

1 

0 


sin 0 

0 

cos 0 


cos 0 

sin 0 

0 ~ 

£3(0) = 

-sin 0 

cos 0 

0 


0 

0 

1 _ 


EULER-ANGLE TRANSFORMATIONS 
\b “ 



cos 6 cos ^ 

sin (p sin d cos \jj - cos 0 sin 0 
cos 0 sin d cos 0 sin 0 sin 0 


cos d sin 0 

sin 0 sin 0 sin 0 + cos 0 cos 0 
cos 0 sin 0 sin 0 - sin 0 cos 0 


-sin 0 
sin 0 cos 6 
cos 0 cos 0 


SOME GENERAL TRANSFORMATION PROPERTIES 


\b- = 

\b' 






^a,c ^c,b 






•^a '•^b 

-la 

*ib 

ia 

O' 

1 

\b = 

ia *-^b 

ia 

*ib 

ia 



i^a-Jb 

l^a 

*ib 

l^a 

'Isb 




Figure 2.2.- Transformations for orthogonal right-handed axes 



’’a.b = E,(»)E 2 ( 0 )E 3 (») 


( 2 . 10 ) 


This is expanded in figure 2.2. The angles determining T^. jj are denoted 

. These, together with suffice to determine all the transformations 

used. Some general properties regariiing inversion and invariance of the transforma- 
tions, and the relation of the axes vectors with the rows and columns of the trans- 
formation are also noted in figure 2.2 for use in the text. These last are obtained 
by noting that the elements of a transformation matrix are direction cosines or the 
dot products among the axes vectors of the two frames. 

The directions of vectors relative to triangle axes are given by the angles 
(e,x). Their definitions and appearance in the triangle axes components of a unit 
vector are shown in figure 2.1; X is the angle out of the triangle plane and e is 
the angle out of the plane perpendicular to the spreader bar (a tilt angle relative 
to the spreader bar). Subscripts indicate the corresponding direction vector; e.g., 
(ei 3 ,X^ 3 >, (c^,X^) correspond to lU, u13, ul- 

In the analysis, it is useful to neglect sufficiently small forces and 
moments. These are considered negligible, provided the corresponding linear or 
angular accelerations are of order 0.01 g or 0.01 rps^ over the operational flight 
envelope; that is, F or M is negligible provided 

I F 1 /m < < 0.1 g 
or 

<< 0.1 rps^ 

where (m, J) are the mass and inertia matrix for the body to which F or M is 
applied. 

Finally, we list the principal assumptions made in the analysis regarding the 
dual-lift system and the flight conditions of interest. 

1 . The level-north reference frame is inertial . 

2. The configuration is arranged as recommended in early studies (refs. 1 

and using a spreader bar and cables (fig. 1.1). Typically, the bridle angle 

5 is in the range 45° to 60° and masses are such that 0 < mj^ << m^^ and m.|,m 2 
are the same order of magnitude as mj^ for the range of loads considered. 

3. The cables are weightless straight lines of fixed length which support only 

tension forces and generate no aerodynamic force. All cables are in tension. 

4. The system is in equilibrium: (a) a reference point on the configuration 

follows the reference trajectory, (b) all redundant variables in the system have 
their commanded values, (c) the load is stable and in steady state attitude, (d) the 
helicopters are in equilibrium. 

5. The maneuver domain of interest is bounded; the angular velocity of the 
path tangent, and the reference trajectory acceleration and Jerk are much smaller 
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than 1 rps, 1 g, 1 g/sec, respectively. The acceleration can be written in terms of 
the path tangent vector and its perpendiculars in the horizontal and vertical 
planes : 


“Y = ♦„!% * TiH 

a = v(cos Yifj - sin cos - Vy(sin yi^^ + cos y^) 

and we assume reasonable upper bounds: 


li 


V cos y|ii»^ 

vi; 

Iv 


< 0.06 rps 

< 0.06 rps 

< 0.15 g 


Heading rate 
Flightpath angle rate 
Speed rate 


< min(0.35 g, 0.06 V cos y) Centrifugal acceleration 


< min(0.05 g» 0.06 V) 
<0.1 g/sec 


Normal acceleration 
Speed Jerk 


^ (V cos yil>y) 


dt 


(v;) 


< 0.1 g/sec 


< 0.05 g/sec 


Centrifugal acceleration rate 


Normal acceleration rate 


6. The angle between the apparent suspended load and the vertical is 

bounded: L (fl,kjj) < S. 

7. The apparent load differences are bounded: 

Ilk - HI, |fl - l£2 - HI < sin 6 

The last two assumptions are conditions derived in sections 6 and 7 which 
restrict flight condition parameters (acceleration and aerodynamic force) suffi- 
ciently to ensure that equilibrium is possible with appropriate geometry; that is, 
with the system right side up and with the helicopters on the side of the spreader 
bar opposite the load with their thrust vectors directed away from the spreader bar 
at any and 6j.. 
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3. FORCE-BALANCE EQUATIONS 


The forces applied at various points on the system are shown in fig- 
ure 3.1(a). The force-balance equations at the load c.g., at the spreader-bar 
endpoints, and at the helicopter c.g.'s are 

F35 u35 + F45 u45 = m. ft 
~ ^ — ■ 

-F13 ull F35 u35 - F34 i^ = -0.5 

-F24 u24 + F45 u45 + F34 i^ = -0.5 f4 (3.1) 

F13 ull - T1 ul = -®i H 

F24 u24 - T2 u2 = f2 

where the apparent loads on the component bodies are 

ft = g + fat - at 
£1 = g + fal - a1 
f2 = g + fa2 - a2 
f3 = fb + ^ 
f4 = fb - Af 
fb = g + fab - ab 

^ it ® (Hb - Mab) 

D 


The force-balance equations for the spreader-bar endpoints are used instead of 
the rigid-body equations for force-and-moment balance about the spreader-bar c.g. 
because the latter do not contain the spreader-bar compression. The development of 
these equations is given in the appendix. The remaining equations are the usual 
rigid-body equations for translational motion of the suspended load and helicopter 
c.g.'s, and they can be written by reference to figure 3.1(a). 

These equations are valid provided no cable collapses (provided F13, F24, F35, 
F45 > 0), and this restriction coincides with the present limits of interest. 

The analysis is carried out principally in triangle axes. The unit vectors in 
equations (3.1) can be expressed in terms of direction angles in these axes as 
defined in figure 2.1. For the present analysis, it is assumed that all the appar- 
ent loads on the right in equations (3.1) can be determined from the reference 
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trajectory. Then equations (3.1), expressed in triangle axes, are solved for the 
following unknown system variables: 

*l’t ®t* '**t Euler angles locating triangle axes 

relative to level-north axes 

F35, F45, F13, F24, F34, T1, T2 Suspension system tensions and 

trust magnitudes 

Tether and thrust direction angles 
relative to the triangle 

Equations (3.1) provide only 15 scalar equations so that 3 of the 18 variables in 
equation (3.2a), or functions of these variables, can be selected independently. 

For these we choose three angles that can be measured and controlled in flight and 
that are later found to be related to factors of engineering interest (see 
fig. 3.1(b)): 

®13 

't 

- *v 

First, for any given flight condition, the thrust-sum, T1 + T2, varies with 
but is very nearly independent of the other two variables. Hence, the tether direc- 
tion can be selected to minimize the required thrust-sum; that is, to preclude 
tether directions which require the two helicopters to tug against each other along 
the spreader bar with mutually canceling thrust components. Second, load sharing 
must be set at a desired value, and this depends on which is the tilt of the 
spreader bar relative to the apparent suspended load, and which is a function of 
♦j., e^, and As will be seen, the bridle-cable load ratio, F35/F45, depends only 
on this angle and, assuming that is selected to minimize the thrust-sum, then 

T1/T2 also depends on only, whereas the sum T1 + T2 is nearly invariant 
with Third, orients the spreader-bar heading relative to the ground 

track. Thrust requirements and load sharing are virtually independent of 0^. In 
manually controlled formation flight |0^| is typically maintained constant in the 
range of 30* to 60*; this value enhances visibility and overshoot safety for the 
trail helicopter which acts to maintain the relative position vector as the lead 
helicopter selects the flight trajectory. In many respects, the situation is simi- 
lar for dual lift, and the same choice of 8^. can be made for manually controlled 
flight. With automatic control, it may be possible to select or vary 0^ over a 
larger range, but a fixed or slowly varying value is expected. A uniquely preferred 
value of 0|- cannot be given here, and results are given for any formation angle. 

Problem Statement (restricted trim): Solve equations (3.1) for the configura- 

tion variables (eq. (3.2a)) in terms of the independent variables (eq. (3.2b)) 


Tilt angle of No. 1 tether from spreader bar 
Spreader-bar tilt angle from apparent load vector 
Spreader-bar heading relative to the reference ground track 


(3.2b) 
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assuming the system masses (m^^, m^, m^, m 2 ) and the apparent loads 

f2jj, f3jj, f4jj) are known, and that the assumptions of section 2 apply. 

This is a restricted trim problem in which the apparent loads are treated as 
parameters in the analysis. This treatment yields closed-form relations for the 
suspension variables and thrust- vector requirements in terms of the apparent loads 
and indicates the general effects of maneuvering and of aerodynamic forces on the 
equilibrium behavior of the system. 

The solution of the restricted trim problem is meaningful provided the apparent 
loads do not depend significantly on the unknowns to be determined from equa- 
tions (3.1); that is, provided each apparent load can be given independently of 
solving the force-balance equation in which it appears. Calculation of the apparent 
loads for any specific system requires numerical solution of the general trim prob- 
lem in which the system's force- and moment-balance equations and constraints are 
satisfied, given the reference trajectory, aerodynamic and engine descriptions of 
the bodies, and assuming equilibrium motion. An algorithm for this purpose is 
sketched in figure 3.2; it is arranged to satisfy the equations, beginning with the 
load and working up to the helicopters, using iteration to account for coupling with 
the rigid bodies higher in the system. Many of the equations required for the 
system below the helicopters in this algorithm are obtained in the present work, but 
its implementation is outside the scope of the present work. 

Figure 3.2 illustrates several points relevant to the restricted trim prob- 
lem. First, the load acceleration, aerodynamic force, and apparent load can be 
computed independently of the remaining system from its rigid-body equations (ignor- 
ing downwash), and are, therefore, functions of the reference trajectory only. In 
the presence of downwash, some coupling with thrust occurs, but it can be assumed to 
be negligible compared with load weight in the restricted trim problem. Second, the 
motion of the triangle is given by kinematic relations from the reference trajectory 
and the selected values of 8^ and after which the spreader-bar c.g. accelera- 
tion, aerodynamic force, and apparent load can be given; they are functions of these 
same variables. Third, the tether-force vectors can be determined from the apparent 
loads of the bodies below the tethers and the selected value of the tether angle, 
e-ig. These forces are then used with helicopter trim routines to calculate the 
corresponding helicopter thrust requirements and apparent loads. Iteration is 
required to account for the helicopter apparent loads in selecting the tether angle 
for minimum thrust. In this step, helicopter aerodynamic specific forces vary with 
thrust, but these variations are weak so that faljj and fa2^j in equilibrium are 
principally functions of the reference trajectory. In addition, these forces are 
small, of the order of 0.1 g or less, so that neglecting their dependence on the 
unknowns of the restricted trim problem will cause little inaccuracy in the 
restricted trim solutions. 

In the next section, the size and nature of the c.g. accelerations and aerody- 
namic forces in equilibrium flight over the operational envelope of interest are 
described before equations (3.1) are solved. 
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4. ACCELERATIONS AND AERODYNAMICS 


The c.g. accelerations and aerodynamic forces of the component bodies of the 
dual-lift system are necessary to the evaluation of the apparent loads in equa- 
tions (3.1). Although their exact calculation from the reference trajectory, the 
aerodynamic descriptions of the bodies, and the assumption of equilibrium motion is 
not attempted in this paper, it is of interest to consider several aspects of these 
quantities before solving the force-balance equations. These are (1) the nature and 
magnitude of differences of the component body accelerations from the reference 
acceleration, (2) the nature and magnitude of their aerodynamic forces, and (3) the 
dependence of these quantities on the unknown suspension variables (eq. (3.2a)). 

The analysis neglects the aerodynamic forces on the cables. Those forces 
acting on the spreader bar are included in the analytical results, but are not 
discussed below. The bar's aerodynamic force is expected to be negligible relative 
to load weights encountered in practice and, therefore, to have little effect on the 
system. It is also remarked that this force can be computed independently of solv- 
ing the spreader-bar force-balance equations (eqs.(3.1)). 

The dual-lift system is assumed to be in equilibrium flight at all points along 
a reference trajectory; that is, (1) a reference point in the system follows the 
reference trajectory, given as (Rojj(t), Vojj(t), aojj(t)) over the flight duration, 
t^ < t < tp; (2) the redundant suspension variables (0^., e^, e-j^) have their com- 
manded values, and their command rates are assumed small (well below 0.1 rps); 

(3) the load attitude and angular motion are stable and in steady state; and (4) the 
helicopters are in equilibrium. 


Accelerations 

In static equilibrium with fixed values of all redundant angles, all parts of 
the system have the saune velocity and acceleration. More generally, this is not the 
case and we examine next the nature and magnitudes of the differences in velocity 
and acceleration from the reference trajectory that can occur in equilibrium flight. 

The position of a point R on a system consisting of rigid bodies and straight 
links of fixed length, such as those in the dual-lift system, can be given as 


R = Ro + ('♦•D 

n 

where ^ is a reference point in the system and (^^”\ n = 1,2,...) is any 
sequence of fixed-length links or line segments connecting ^ and R. The inertial 
velocity and acceleration of R can be given from Coriolis' equation specialized to 
this case: 
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(4.2) 


n 


a = 


where is the 

equation (4.2), the 


6V = |V - Vo| 
6a = la - ao| 


L r'(n) .o(n) (n) , (n) *d('^)m 

[w ® AR **■ H ® ® — 'J 

n 

angular velocity of relative to inertial space, 

differences from the reference trajectory are 

= < max{|( 0 ^”^ I (sin Z.(ai^”\ ^^”^)(}L 

n n 

iV^/’’(n) AD(n) (n) (n) ^ AD(n)\i 

n 


From 


(4.3) 


n 

< max{|w^"^ + 0 )^”^ (sin AR^^^) | }l 

n 


where L is the maximum dimension of the system (about two rotor diameters). These 
differences would be zero if the system were a point mass (L = 0) or if there were 
no angular motion = 0, n = 1,2,...), and they would be negligible for suffi- 

ciently small angular motion and system dimension (that is, if 

2 

max{|u^^^| + }L << 0.1 g) 

n 

We omit a detailed application of equations (4.2) to the dual-lift system, but 
the following are noted. The angular motion of the system is determined from the 
rates of the independent angles in the problem acting in the directions about which 
each rotation occurs. These are the direction rates of the reference velocity 
(y, i ) and apparent suspended load, the command rates of the redundant suspension 
variables (B^, e^o)> ^"d the load and helicopter attitude rates for equilibrium 

flight (i, a, B for each body). On quasi-steady flight segments (constant speed or 
accelerating straight line and steady turn segments) all of these angles except 
i|>^ can be assumed invariant or to have very small rates, so that the angular motion 
of all parts of the dual-lift system is either negligible or dominated by the veloc- 
ity heading rate in turns; 

ii ~ *^v^ ’ n=1,2, ... 
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For these flight segments » the system approximates a rigid body with uniform angular 
motion and the velocity and accelerations of its parts differ from the reference 

• *2 

values, Vo, ao, to terms of order I’^ylmax^' respectively, where L is 

the maximum system dimension and U I is below 0.1 rps for the trajectory 

* v'max 

domain of interest previously stated. During brief transient maneuvering between 
steady-flight segments (turn entries and exits, v and y - transients), the direction 
rates of the apparent loads contribute to the angular motion of the system, the 
system's bodies change their relative positions, and the maximum angular velocity 
can be greater, but is still expected to be small, of the order of 0.1 rps for the 
acceleration and Jerk limits previously stated, with corresponding maximum accelera- 
tion differences in equations (4.3) under 0.1 g. 


Helicopter Aerodynamic Forces 

The nature and magnitude of the helicopter aerodynamic forces exclusive of the 
rotor forces can be evaluated approximately from trim solutions for the isolated 
helicopter. Aerodynamic models based on available wind-tunnel data (e.g., 
refs. 12-14 for the HO-4 and UH-1H helicopters; ref. 15 for the CH-54; ref. 16 for 
the UH-60; and ref. 17 for the CH-47B) along with trim-solution algorithms for 
isolated or externally loaded helicopters in static equilibrium flight (e.g., 
refs. 12,16,18) are available from simulation programs. The trim algorithm asso- 
ciated with the simulation described in references 13 and 16 was used to obtain data 
for the UH-60. 

Contour plots of the aerodynamic specific force in level-heading axes aligned 
with the heading of the air-velocity vector are given in figure 4.1. Results are 
given for two weights corresponding to a helicopter weight of 12,000 lb with both 
light and heavy loads (2,000 and 10,000 lb) attached at the helicopter e.g. Several 
features can be seen. First, the load weight and thrust have a minor effect and it 
is principally on the vertical component. Second, the magnitude of the longitudinal 
component is small and under 0.1 g over the speed range shown (at low speeds it is 
forward (positive), and otherwise increases principally with Ve^), and the lateral 
component is virtually negligible everywhere (a discontinuity occurs in the contours 
at Ve = 60 knots where the trim algorithm switches from selecting the helicopter's 
redundant variable as S = 0 to <|) = 0). The vertical component is also small, 
with similar maximum magnitude; it provides little aerodynamic lift (at higher 
speeds it varies principally with y and is downward (positive) at higher values 
of y)* And third, the plots of the specific force magnitude |fa| confirm that 
the helicopter's total aerodynamic specific force is generally small (below 0.1 g 
over most of the domain shown) and that it can, therefore, be expected to have 
little effect on the equilibrium rotor-thrust requirements or other system 
characteristics . 

These data are specific to the UH-60 helicopter in static equilibrium with the 
load attached at the e.g. However, additional calculations for the tandem-rotor 
CH-47B with and without a 15,000-lb slung load, using the simulation described in 
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Figure 4.1.- UH-60 aerodynamic specific force in static equilibrium flight, 

Fau/W^,: Wj, = 12,000 lb. 






reference 17, gave the same trends, and those trends are expected to characterize 
other helicopters and load attachment geometries generally. Thus, the helicopter 
aerodynamic specific force depends weakly on the thrust magnitudes and tether 
forces, but this dependence can be neglected in solving equations (3.1) since these 
forces are generally small (of the order of 0.1 g) over the velocity range of inter- 
est; as a result. Inaccuracies in the present analytical results owing to the 
neglect of this dependence are not significant. Further, the helicopter aerodynamic 
forces have no effect on the system below the tethers in the results derived later. 


Load Aerodynamic Forces 

The physical characteristics of external loads that might be transported by 
helicopters vary substantially and the loads are of uncertain stability in forward 
flight before system testing. In general, a given load and suspension system will 
be stable beneath helicopters in static equilibrium flight over some speed range. 
This stability-limited speed range can be determined empirically in wind-tunnel and 
flight tests, and it can be significantly smaller than the power-limited speed range 
of the helicopters depending on suspension-system design properties (restraints in 
yaw, pitch, and pendulum motion, and characteristic lengths) and load properties 
(aerodynamic force and moment characteristics, weight, c.g. location). Test data 
are available for single helicopter suspensions (refs. 8-10,19-21, and the 1976 
survey article, reference 22), but not for the dual-lift suspension of this study. 
Here, the lack of yaw restraint and long pendulum length are unfavorable factors, 
and the lower bound on load weight (assumed to exceed the single-helicopter payload) 
is favorable for a given load shape and size. It is beyond our purpose to evaluate 
the load-speed range combinations that satisfy the load stability assumptions of 
this paper. However, it is reasonable to assume for the present that the dual-lift 
and single-cable suspensions have similar load stability boundaries, since both lack 
yaw restraint. This boundary is noted in references 10 and 19 as being 40 to 
65 knots for various difficult but typical loads (e.g., boxes, vehicles, and bridge 
sections); reference 21 identifies four general categories of military loads and 
notes maximum operational speeds of 40 to 100 knots, depending on category. 

Measurements of the static aerodynamics of some typical loads are given in 
references 23 and 24 (containers of sizes 8x8x8 ft, 8x8x20 ft (MILVAN), and 
8x8x40 ft; a bulldozer; and the M109 self-propelled 155-mm howitzer). These loads 
are convenient examples for the examination of maximum aerodynamic specific forces 
and exemplify trim behavior for a large group of loads. For brevity, we omit a load 
trim analysis based on these data and infer trim attitudes and maximum stable speeds 
from the literature on single helicopters with single cable suspensions (ref. 22). 
The elongated boxes of reference 23, and other elongated loads with a definite major 
axis and weak directional stability when aligned with the wind, generally stabilize 
broadside to the free stream in a maximum drag heading with small or negligible lift 
and side forces. The cubic box, and other loads which are axlsymmetric about the 
cable, generally rotate continually over their stable speed range. The aerodynamic 
force is again principally drag, which is nearly independent of yaw angle. Last, 
the M109 is a dense load with some degree of weathervane stability at the minimum 
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drag heading. Flight observations are not available, but a similar load (the Ml 10 
self-propelled 203-mm howitzer) trims longitudinally up to 120 knots, and then trims 
broadside for the remainder of its stable speed range. 

The wind-tunnel measurements in references 23 and 24 are given for force/Q, 
which is assumed independent of Q. The equilibrium aerodynamic specific force on 
the load can be written as 


fat 



W,ia 


(4.4a) 


where fat and i^ are the magnitude and direction of the static aerodynamic force 
vector, and vary with Oj^, The trim force is dominated by drag for all the 

example loads and other loads with the same trim attitude behavior: 


fat 


-f-'l - 1 


(4.4b) 


where D, are the load’s drag and air velocity direction, respectively. The 

loads can be characterized by the proportionality factor, ( 1/W)(D/Q)^.^^j„, which is 
the drag specific force per unit dynamic pressure, and indicates how rapidly spe- 
cific force increases with airspeed. Data for the example loads are given in 
figure 4.2, including maximum stable speeds (critical speed), operational speeds, 
and the corresponding specific forces. These speeds are typical or approximate 
values from wind-tunnel or flight tests for the same or a similar load. 

The MILVAN’s maximum stable speed appears to be independent of weight, so that 
the maximum specific force to which it is subjected decreases with weight along a 
straight line (fig. 4.2(b)) over the range 0.145 g to 0.525 g. Maximum specific 
forces for the other loads fall within this range. The Ml 09 has a much lower drag- 
specific-force rate than the MILVAN and is subjected to much lower specific force at 
a given speed, with values below 0.1 g at all speeds below 100 knots, and corre- 
spondingly less influence of load aerodynamics on the system equilibrium. Specific 
forces at operational speeds are also shown. Since the same speed is assigned to 
all containers, these fall along a straight line in figure 4.2(b) and are below 
0.3 g for the MILVAN at all weights. 

These data give the maximum specific forces for the example loads cited above 
and indicate bounds for other other similar loads with drag-specific-force rate in 
the same range. The discussion also indicates that most stable slung loads are 
subjected to aerodynamic forces that are dominated by drag and that are, at moat, a 
moderate fraction of 1 g. A load is unlikely to be stable if the specific force 
approaches 1 g. Further, later sections of this report indicate that unacceptable 
equilibrium geometry is required if the combined load aerodynamic force and acceler- 
ation is sufficiently large (approaches 2 g). 
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/D\ 

'Q'TRIM' 

ft2 

W, 

lb 

' 

'^WtrIM' 

(Ib/ft2)-1 

^CRITICAL' 

knots 

'fa 'max' 
9 

VqP'L' 

knots 

'fa'oP'L' 

9 

M109 

85 

43,500 

0.003 

140 

0.20 

75 

0.06 

8X8X8' BOX 

79 

2,500 

0.032 

65 

0.46 

45 

0.22 



9,000 

0.009 

90 

0.25 

45 

0.06 

8 X 8 X 20' BOX 








EMPTY 

216 

5,000 

0.043 

60 

0.525 

45 

0.30 

FULL 


21,000 

0.012 

60 

0.145 

45 

0.08 


(a) DATA FOR EXAMPLE LOADS 



(b) EQUILIBRIUM SPECIFIC FORCE 

Figure 4.2.- Maximum equilibrium load aerodynamic specific force (refs. 22-24). 
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5. RESULTS FROM THE LOAD EQUATION 


The following results are obtained from the first of equations {3.1)t 


^.fi = 0 

F35 = 0.5 cos(6 - ej^)/cos 6 sin 6 

F45 = 0.5 cos(i5 + t^)/cos 6 sin 6 

P|j = F35/F45 = cos(6 - ej^)/cos(« + 
F35, F45 >0 IFF |ej^| < n/2 - 6 
F35 = F45 IFF e = 0 

jL 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


These results readily follow from analysis of the dot products of the load equation 
with i^, after expressing the cable directions, u35 . u45 . in triangle axes 
(refer to fig. 2.1). These results show the following. 


First, since is perpendicular to the triangle it follows from equa- 
tion (5.1) that the specific apparent suspended load ft is necessarily in the 
plane of the triangle at all times. In that case. 


X = 0 

and (5.7) 

ut^ = (-sin 0, cos e^) 

This reflects the assumption that the cables are straight lines which support no 
side force. As a result, is the pitch angle of the spreader bar from the plane 
perpendicular to fa; the spreader bar is level relative to fa when = 0. 

Second, it follows from equations (5.2) and (5.3) that 


F35 > 0 IFF -ir/2 + 6 < e < ir/2 + 6 


F45 > 0 IFF -ir/2 - 6 < e < ir/2 - 6 

Xr 


Combining these gives the result in equation (5.5) for the range of spreader-bar 
tilt relative to fa for which the bridle cables do not collapse. For 6 = 60®, 
this is a range of ±30*. At the extremes of one of the bridle cables is paral- 

lel to fa and carries the entire suspended load, mj^fa, whereas the other col- 
lapses. It will be seen later that the tether cables do not collapse at any so 

that equation (5.5) gives the range of for which equations (3.1) are valid. 

Third, the bridle-cable tension ratio, equation (5.4), depends only on the 
spreader-bar tilt, Equal tension requires that the spreader bar be 
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perpendicular to the apparent suspended load (be level relative to fa). It will be 
shown later that, assuming optimum choice of tether angle, the tether- tens ion ratio 
is very nearly the same as Pj^, and equal loading of the helicopters requires that 
be nearly zero. 

Last, the variations of cable forces and force ratio with e,,6 are shown in 
figures 5.1 and 5.2. 



Figure 5.1.- Pendant cable force ratio: p = F35/F45. 

P 



Figure 5.2.- Pendant cable forces. 
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6. BRIDLE TRIANGLE ATTITUDE 


The triangle attitude is defined by its body-axis Euler angles, 
where 9^ are the spreader-bar heading and pitch, and <|>j. is the load roll or 

swing angle about the spreader bar from the vertical plane. These angles can be 
given in terms of >1)^, 8^, and uHj^ (the direction of ft) from the following 

identities : 


and 


or 


(•t = 




( 6 . 1 ) 


( -sin 
cos 







( 6 . 2 ) 


where (ut^j, uft , ut^) are the known components of ut^. Equation (6.1) gives i|»(. 
in terms of il»^, 6f Equation (6.2) imposes the conditions that ut be in the 
triangle plane and that it be tilted, relative to the spreader bar, by the angle 
Gj^. This yields the following relations for the spreader-bar pitch and load swing 
angles (derivation omitted): 


sin (t)^ = (ua^ sin 8^. - cos 8^.)/co3 
sin (0^ - E) = sin cos 


(6.3) 

(6.4) 


where 


tan z = (ul^ cos 8^ + sin 

Equations (6.1), (6.3), and (6.4) give the load swing and inertial spreader-bar 
attitude to be maintained by the helicopters for equilibrium flight as a function of 
flight condition (ulf^, tp^) and the selected values of 8^.. 


Solution Existence Conditions 

Equations (6.3) and (6.4) yield a solution for (p^., 0^., provided their right- 
hand sides have magnitudes less than 1. Necessary and sufficient conditions for 
this are 

2 2 1/2 

((p^,0^) exist IFF [uJ,^ + cos 8j. + ul^ sin 8^.) 1 > |sin e^| (6.5a) 
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Further restrictions are needed to ensure that the triangle is right side up 
(|(>tl» |6j.| < tt/ 2) at any choice of 6|.: 

exist at any 6^. and ^ ^ (6.5b) 

or at any choice of e^f 

(♦j.»0^.) exist at any and le^^l < ir/2 - 6 and |<|)^|,|9^| < it/2 IFF ut^ > cos 6 

(6.5c) 


These are progressively more restrictive conditions on ut which allow for 
satisfactory solutions over progressively larger domains of the independent varia- 
bles. It is expected that the narrowest condition, equation (6.5c), is satisfied in 
dual-lift operations; that is, that ut^ > cos 5. Then a 
solution exists for any choice of 6|. and any admissible 
with the triangle right side up. Note that is defined by 

the angle between the apparent suspended load and the true 
vertical, denoted in the sketch, and that the condition in 

equation (6.5c) can be written as 

or 



'-H 




fl-kN 


< tan 6 


The angle depends entirely on the combination fan - at: it increases, for 

example, with load drag, centrifugal acceleration, or load lift (provided ft • ijj 
or fn • is nonzero), and it reaches 60® if, for example, the horizontal compo- 
nent of the apparent load is 1.73 g. This condition provides general limits on 
an and on the load aerodynamic specific force for which it can be appropriate to 
use the present mode of transport. It is expected that within the stated limits of 
interest in an, and for load aerodynamic specific forces encountered with stable 
slung loads (principally moderate drag force with little lift) that the condition in 
equation (6.5c) is met, and this is assumed to be the case hereinafter. 


Results 

To examine the variation of the triangle pitch and roll attitude with flight 
condition along a trajectory, it is convenient to express the direction un^j, in 
terms of direction angles. This can be done in several ways. 
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First, can be given in terms of where 



/sin 003 8„A 

3ln 3ln 6„, 

\ 5 , / 


( 6 . 6 ) 


or, inversely, 


sin 5^ = (ua^ + 


0 < < it/2 


2 2 1/2 

(sin 0^^, cos = (ua^,uay)/(ua^ + ua^) 


■TT < < IT 


Here, is the angular offset of ua from the vertical, as sketched above, and 
indicates the net effect of trajectory acceleration or significant aerodynamic force 
on the direction of the apparent suspended load. It can reach magnitudes of 20 ® to 
30 * for the regime of interest here. The angle 0^^ is the heading of fa and 
faa - ^ relative to the ground track. Although 0^^ is undefined when fa is 
vertical (when = 0) there is no computational difficulty because are 

independent of 0^j^ in this case. The results in equations (6.3) and (6.4) can now 
be written as 


sin 4)^. = sin sin(0^ - 0^j^)/cos 

sin(0^ - 2) = sin cos r/cos 

tan Z = tan cos(0^ - 0^^^^) 


< ff/2 
(e^ - z| < n/2 

III < tt/2 


(6.7) 


These give the triangle attitude in terms of three variables, c 


suspended load. 


®t - 


is the spreader-bar heading relative to the heading of tne apparent 




and 


The required spreader-bar pitch, equation (6.7), has the following properties: 


|i| ^ h 

3,.r = ° 

(e + A ifej^^tO 


where 

|sin A| = ]sin cos Z/cos > jsin 

The required load roll about the spreader bar in equation (6.7) can be made 
zero at all times along a trajectory by continually aligning the spreader-bar 
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heading with that of the apparent load and varying the spreader-bar pitch with 
tal - at as required; that is, 

h = ♦t = ° \ = h 

Similarly, in the case of equal bridle-cable loading, e^. can be made zero at all 
times by aligning the spreader bar perpendicular to 8^^; and then the load swing 
varies with faft - as required. That is, 

if: = 0 and 8^. = ± ir/2 , then: ~ ° 

Thus, it is possible to fly the system such that either <j»(. or 0^. is invariant 
along a trajectory by maintaining a fixed spreader-bar heading relative to Byj^. 
However, note that 8yj^ undergoes large changes at turn entries and exits and at 

other transitions between steady flight segments so that flying the system with 

invariant load swing or spreader bar pitch appears impractical. Plots of the map 
e ) (<l>^,0^) are given in figure 6.1 for = 0 (pj^ = 1) and 

= '0* (pjj = 1.9 if 6 = 60®). The trends noted above are evident for the lines 

corresponding to |8. - 8 *|e {O, ±90°, 180°}. 

Second, an alternative form for ua^j can be given in 
terms of pitch and roll angles relative to level-heading 
axes, (see sketch) where 0^ is a rotation 

about Jjj from the vertical which reflects acceleration 
and aerodynamic force along the ground track and is a 

roll angle which reflects acceleration and aerodynamic force 
lateral to the ground track. For these angles: 



or, inversely: 



/cos sin 9^ \ 

-sin 1 

V cos <t> . cos 0 . / 
\ uft ul' 


sin 

d» . = -uA 
^ut y 


< ir/2 

tan 

0 . = ut /uS, 
Ut X z 

IVI 

< ir/2 


( 6 . 8 ) 


The range of these angles encountered in quasi-steady flight over the operational 
domain is examined in figure 6.2(a) for the MILVAN, with equilibrium aerodynamic 
force calculated from equation (4.5). First, results are shown for speeds and turn 
rates over their ranges of interest for the MILVAN, with the remaining flight condi- 
tion variables (y» v, wind) taken as zero. The angle 0^^ decreases from 0® to 
-10® with airspeed because of load drag, and varies with turn rate over a 

range that increases to ±14® with speed. Second, the remaining variables expand the 
ranges of encountered in quasi-steady flight to the boundary shown in 

figure 6.2(a) when taken over their maximum ranges of interest. Flightpath angle 
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0^, deg 



0t, deg 



Figure 6.1.- Triangle attitude: vs 


30 








Figure 6. 


BOUNDARY FOR QUASI-STEADY FLIGHT 



(a) DIRECTION ANGLES OF APPARENT LOAD 



Apparent-load direction angles and triangle attitude in quasi-steady 
flight for the MILVAN = 15,000 lb). 
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has no significant effect, whereas v moves the figure left or right, and the wind 
acts to change the drag vector, depending on wind direction and magnitude; these 
effects collectively expand the direction angle ranges of interest to -20® to 10° 
and ±20° for and 41^^, respectively. Last, it is noted that the apparent-load 
magnitude, ft (not shown), is insensitive in these variations in flight condition 
and is within the range of 0.99-1.1 g for the entire set of flight conditions. 

The triangle attitude angles can be given from equations (6.8), ( 6 . 3 ), 
and (6.4) as 


sin = (cos (|> sin 0 sin 8,. + sin 4) cos s,.)/cos e 

C UJt U3C u UjL w 

tan 2 = tan 0^^^ cos 6^. - tan 4>^^^ sin e^/cos 0^^^ (6.9) 

sin(0^ - Z) = sin cos E/cos 4*^^ cos 0^^^^ 

This form of the results is convenient for examining the variation in triangle 
attitude along a trajectory when the system is flown with fixed formation angle, 

6^. Plots of the map given in figure 6.3. For the 

longitudinal formation (8^. = 0), the map is an identity in the case of equal bridle- 
cable loading (ejj = 0); the load roll angle changes with turn rate and aerodynamic 
force lateral to the ground track, while the spreader-bar pitch is nearly fixed, and 
the spreader-bar pitch changes with speed and longitudinal aerodynamic force while 
the load-roll angle is nearly fixed. The reverse trends occur for the transverse 
formation; that is, spreader-bar pitch changes for turns (outside helicopter up) or 
side force, and the load swings backward/forward to accelerate/decelerate or in 
response to load drag. For the 45° formation angle, spreader-bar pitch and load- 
roll angle both change for turns and speed changes or as a result of load drag. 

This is further illustrated in figure 6.2(b) which shows the triangle attitude for 
steady turns and the boundary of the region required for quasi-steady flight in the 
case of the MILVAN with 8^ = 45° and = 0. This boundary is defined as shown 
in figure 6.2(a) and encloses a range of -26° to 8° in both 0^j^ and 4>u)f other 
values of 8^ and the figure rotates and translates on the plot with correspond- 
ing changes in the angle range required in quasi-steady flight. 

All of these variations of triangle attitude with maneuvering and aerodynamic 
forces are defined by the basic properties stated in equations (6.1) and (6.2); 

(1) the triangle necessarily contains f^, (2) the spreader-bar tilt relative to fa 
is 6j^, and (3) the spreader-bar heading from the ground track is 8^.. 
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Triangle attitude map: {<l»4.,64.) vs (♦ ,0 
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7. SPREADER-BAR AND HELICOPTER EQUATIONS: GENERAL RESULTS 


The tether forces can be obtained from the spreader-bar equations in equa- 
tions (3.1) and from the previous results for F35, F45 (eq. (5.1)). First, it is 
convenient to separate the bridle-cable forces into a term parallel to fa and a 
remainder: 


F35 u35 = y m^^fa - A 

F45 u45 = (1 - v)nij^fl + A i^ 


(7.1) 


where 


y = 0.5 COs(6 - 6j^)/COS 6 cos 

A = 0.5 m.fa cos(6 - e,)cos(a + e )/sin 6 cos 6 cos e 

Sr Xr Sr 

The terms with coefficients y, 1 - y, are the apparent suspended load, fa, distrib- 
uted to the endpoints of the spreader bar in the same ratio as the bridle-cable 
tensions; that is, 0 < y < 1 for in its admissible range, and 

y/(1 - y) = p|j. The remaining terms in A are forces along the spreader bar in 
opposite directions, and produce the spreader-bar compression caused by the bridle- 
cable forces. Using equations (7.1) in the spreader-bar equations yields the tether 
forces: 


F13 u13 = F13° - B 
F24 u24 = F240 + B i^ 


(7.2) 


where 

F13 ° = 0.5 mj^f3 + y m^^fa 
F24 ° = 0.5 mj^f4 + (1 - y)m^fa 
B = F34 + A 

Here, F13 ° . F24 ° are the tether forces for the special case that the spreader-bar 
compression is due entirely to the suspended load (by means of the bridle cables) 
with no effect from the helicopters (by means of the tethers); that is, for the case 
in which B = 0 and F34 = -A. The superscript ( )® will hereinafter refer to 
quantities corresponding to this special configuration, which can be obtained by 
appropriate choice of tether angle. 

Equations (7.2) provide six scalar equations containing seven unknowns (F34 
or B and the tether-force components). The analysis is simpler if B is treated 
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as the Independent variable, but this is equivalent to treating F34 or 
as independent. 

The tethers can always be placed in tension at any given B by locating the 
helicopters relative to the spreader bar such that 

u13 = (FI 30 - B i^)/|01° - B ij-l 

u24 = (F240 + B i^)/lF24° + B i^.| 

and then FI 3, F24 > 0. However, interest is limited to the case in which this can 
be done with the helicopters on the side of the spreader bar opposite the suspended 
load; that is, to cases for which uH*k. , u24*k^ > 0, or, equivalently, 

F13 °-k^.. F24 °«k^ > 0. From equation (7.1), 

if; f3j^, f4j^ > 0 , then; F13^, F24° > 0 for all |e^| < f - « (7.3) 

Attention is called to the notation ( )^, ( )j, ( )j^ to denote the t-axes compo- 

nents of a vector here and in the remaining analyses; for example, 
f3(. = (f3j^,f3j ,f3jj) = (£3'it’^‘-t»^^’-t^ • ^ simpler, more conservative condition 

than equation‘s (7. 3) can be given by writing f3jj,f4|^ in terras of the apparent-load 
differences; 


f3jj = (f3 - fi + ft ut) • k^ = (n - H) • Kt 

f4j^ = (f4 - fa + fa ua) • k^ = (f4 - ft) • cos 

and then 

lO - lil Ifi “ lil 

if; < sin 6 , then; F13g,F24° > 0 at all |ej^| < f - 5 

(7.4) 


Or even more conservatively. 


if; 


|ab - at| + 


Ml + Ifat - fab| 
“Ta 


< sin 6 


\ 


then; F13g,F24^ >0 at all hjil < f - « / 


(7.5) 


In the present context |ab - aa|> iMl snd if the aerodynamic force differ- 

ence satisfies | faa - fab I < fa sin 6 - |ab - aa| * (that is, is below about 
0.7 g) then the conditions in equations (7.3) to (7.5) are met. In that case, 
equation (7.2) yields the following results for the tether angles relative to the 
triangle and the tensions; 
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tan = (B - F139)/(F13°^ + 

P P 1 /P 

tan e^n = -^B + F24p/(F240‘^ + F24° ) ' 


tan X ^3 = F13o/F13g 

tan = F240/F240 

2 2 21/2 
F13 = [(F13? - B)*^ + F13J + F13g ] 

F24 = [(F240 + B)^ + F24j^ + F24o^]^^^ 


(7.6) 


These solutions for the angles follow from their relations with components of 
u13^»u24^., defined previously in figure 2.1. All angles are below 90", and 
F13,F24 > 0, provided the conditions in equations (7.3), (7.4), or (7.5) are met. 
The above results are listed in the equation summary, figure 7.1, along with alter- 
native expressions for some variables in terms of ei3»£24* Either of these angles 
or F34 or B can be taken as the independent variable. 


The required helicopter thrust vectors in triangle axes can now be obtained 
from the last two of equations (3.1) and the above results for tether forces: 


where 


X1 = II® + B i^ 

► 

T2 = T20 - B i^ 

11° = -F13° - m^n 

T20 = -F240 - m2f2 


The thrust-direction vectors ul, are defined as 

ul = -n/T1 ' 
u2 = -T2/T2 


(7.7) 


(7.8) 


The negative sign is adopted in equations (7.8) so that the direction angles of ul 
u2 relative to triangle axes will tend to be small. To obtain angles below 90" 
(thrust directed away from the spreader bar) we require T1jj,T2j^ < 0, and this is 
obtained at all admissible provided f1[^,f2|j,f3(f ,f4j^ > 0, or, more simply and 

more conservatively. 
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DEFINITIONS: 

H = 0.5 cos(5 - ej)/cos 5 cos eg 

F13° = 0.5 mjj f3 + ju mg f£ 

F24° = 0.5 m|3 f4 + (1 - m) mg fg 

T1° = -F13°-m.| f1 

T2° = -FM° - m2 f2 

2 2 

ai = (TlP +Tlg ) 

2 2’/* 

* (T2P + T2g ) 

2 2 

(F13P + F13g ) 

2 2’^ 

a4*(F24P +F24g ) 

A = 0.5 mg f£ cos(6 - eg) cos(6 + eg)/sin 8 cos 5 cos eg 
B = F34 + A 

ASSUMPTIONS: 

f2^, f3^, f4^, F34 OR e,3 OR e24« mg, mj,, m^, m2, 8 ARE KNOWN 

|eg|<7r/2-6 

fbk, f Ik, f2k > 0 OR (MORE CONSERVATIVELY) If1 - fgi, |f2 - f£|, |f3 - fg|, |f4 - f£| < « sin 8 

TETHER DIRECTIONS AND FORCES: 

F13u13= F13°- Bi^ 

F24 uM = F24 + Bi^ 

1^13!/ 1^24^' 1^24^ 
tanei3= (B- F13P)/«3 

tan €24 “ + F24P)/tt4 = -{F13P + F24f + ^3 tan ei3)/a4 

tanXi3» F13P/F13g 
tan X24 = F24P/F24g 


Figure 7.1.- Equation summary; spreader bar, tether, and thrust forces. 



2 n ^ 

F13 = [(B - F13f) + tt3^] = a3/cos 6^3 

2 , 

F24 = [(B + F24P) + 04^] = a4/cos €34 

SPREADER-BAR FORCE: 

F34 = -A + B 

F34 = «3 tan 613 ~ 0.5 ni£ f£ cos(6 - ej)/sin 5 + 0.5 mj, f3j 
F34 = -^4 tan €34 - 0.5 mg f£ cos(5 + eg)/sin 6 - 0.5 mj, f4j 

THRUST DIRECTIONS AND MAGNITUDES: 

XI = -0.5 mjj f3 - mg f£ - m., f1 + Bi, = T1° + Bi^ 

12 = -0.5 mj, f4 - (1 - m) mg f£ - m3 f2 - Bi^ = 12° - Bi^ 

X2 + T2 = -mg f£-mjjfb-m^f 1 -m 3f2 

u1 = -I1/T1 

u2 = -T2/T2 

kil» l^ilr 1^2^' ^^2 ^ *^Jt /2 

tan e.j = (T1? + B)/a.| = (a3 tan €^3 - m.j 

tan 63 = (T2° - B)/o(3 = («4 tan €34 - m3 f2j/a3 

tan X-, = TlP/Tlg 

tan X3 = T2P/T2g 

2 o 

T1 = ((Tif -t- B) +a,^] =ai/cosei 
2 

12 = [(T2P - B) + 03^] = «3/cos €3 

2T = T1 +T2> Img ffi + m|,fb + mi n + m3f2| 

, T1 _ «lCOse2 

Ot = “ 

' T2 «3 cos 

F igure 7.1.- Concluded . 
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(7.9) 


if: 


III - Hi IH - Hi IH - fil m - fa| 


fa » fa ’fa ’ fa 

then: T1. ,T2. <0 at all |e.| < x - 6 


< sin 6 , 


These results, along with expressions for the thrust direction angles and magni- 
tudes, the thrust ratio, and thrust-sum are listed in figure 7.1. Alternative 
expressions in terms of the thrust angles given for some variables, 

and these show many similarities of form to the tether force results. 


The following remarks pertain to figure 7.1. 

1. The out-of-plane tether angles, 
fa, fb) and on the value of and are 

nij^ f3, f4 << ro^^fa then ^i3»^24 
possibly near the extremes of where 

tan ^3j/f3jj 

In general, the tethers lie in the plane 
windless hover) and otherwise lie out of 

2. The spreader-bar force is given 

direction angles, ^24’ conversely. Any of the three can be regarded as 

the independent variable. The expression for F34 contains three terms owing to, 
respectively, the tether direction, the suspended load, and the spreader bar. At 
the extremes of ej^ (at |ej^| = it/2 - 6), F34 is independent of the load. 

3. The expressions for the required thrust, T1^,T2, have four terms which show 

that each helicopter carries its own apparent load plus very nearly half the 
spreader-bar apparent load (f3, f4 differ negligibly from fb), plus a fraction of 
the apparent suspended load determined by plus mutually canceling terms of 
arbitrary magnitude along the spreader bar. At the extremes of one helicopter 

carries the entire suspended load, and the terms in for the other are 

zero. At = 0, the suspended load terms are the same for both helicopters. 

4. The out-of-plane thrust angles, X^,\ 2 , are independent of B or F34 or 

tether direction angles. In general, the thrust vectors lie in the triangle plane 
if + 2m^f1j = ffl|3f4j + 2 m 2 f 2 j = 0 (e.g., at windless hover) and otherwise lie 

out of the plane, principally to balance fij or f2j. 

5. The thrust sum and corresponding fuel requirements vary with B and can be 
minimized with respect to B. Solutions for this optimum configuration are given in 
the next section. At other values of B, only the mutually canceling components 

of T1^,T2 along the spreader bar (along i^.) vary. 

6. The vector sum T1_ + ^ is invariant with B and equals the sum of all 
apparent loads, as is expected from the force balance of the entire system. There- 
fore, the magnitude of the apparent-load sum is a lower bound for the 


^13’ ^24’ trajectory (on 

independent of F34 or B. Since 
small angles at all le^^l < ir/2 - fi, except 
one of the tethers is nearly unloaded, and 

tan X, 


or 


'24 




of the triangle if f3j = f4j = 0 (e.g., at 
the plane by small angles. 


in terms of B or in terms of the tether 
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thrust-magnitude sum, zT. The lower bound is obtained only if the system can be 
flown with Tl^ and T2 parallel to the apparent-load sum. 
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8. SPREADER-BAR AND HELICOPTER EQUATIONS: TETHER-ANGLE OPTIMIZATION 


The tether angle equivalently, 624 or F34 or B) is an independent 

variable in the general results (fig. 7.1) and can be selected to optimize some 
dependent function, such as thrust margin or fuel rate. The system thrust margin 

Mr = - T1) + (T2^^^ - T2) = (T1 + T2^^^) - Z„ (8.1 

X nicLX m^x nidx x 

is maximized at the minimum thrust sum. 

Fuel-rate dependence on thrust was examined empirically for the UH-60 as an 
example case. First, test-stand data for its engines (two GE-T700s) indicate that 
fuel rate is very nearly linear and increasing with power output, and is accurately 
modeled in the form 


w^ = a + bP 

where a and b vary with ambient air conditions and Mach number (see fig. 8.1). 
Second, the required power is nonlinear in thrust, but the nonlinearity is usually 
weak over the thrust range of interest. A simulation model of the UH-60 (refs. 13 
and 16) was trimmed in various static equilibrium flight conditions for a range of 
weights in order to obtain power required versus thrust output. The results 
(fig. 8.2) show that the required power increases with thrust at all flight condi- 
tions except in the region of high-speed descending flight where it is insensitive 
to thrust. In addition, it is adequately approximated as linear in thrust over the 


Wf = a + bP 

a/(«\/F)0® = 146.1 Ib/hr 

b (1 + 0 . 2 m 2 ) 1-6 = 0.360 Ib/hp-hr 

5 = PRESSURE RATIO (AMBIENT/2116 Ib/ft^) 

d = TEMPERATURE RATIO (AMBIENT/518.688 °R) 



Figure 8.1.- Fuel rate versus power: GE-T700 engine. 
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THRUST, lb THRUST, lb 

Figure 8.2.- Power required versus thrust in unaccelerated flight: UH-60. 


range of flight conditions and thrust shown in the figure; that is, the required 
power for the UH-60 is of the form 


P = c + dT 

where the coefficients c and d vary with flight conditions. Then, fuel rate and 
thrust are linearly related: 

Wj. = A + BT 

The same degree of linearity and trends for the slope d are expected to character- 
ize other engines and helicopters so that, in general, the dual-lift fuel rate sum 
is approximately linearly related to thrust: 

iw^ = (A^ + A2) + B^TI + B2T2 ( 8 . 2 ) 
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For flight conditions with positive and B 2 , the minimum fuel is obtained by 
minimi 2 ing the thrust-sum (weighted thrust-sum) in the case of identical (nonidenti- 
cal) helicopters. 

We select the tether angle to minimize the thrust-sum in this study. This 
choice yields the maximum system thrust margin in all cases and, for Identical 
helicopters, it is expected that the corresponding fuel requirements will be mini- 
mized or insensitive to thrust and tether angle at all flight conditions. In addi- 
tion, results are given for the configuration in which the apparent suspended load 
is distributed to the two helicopters in the same ratio as the bridle-cable force 
ratio; this yields simple general results for tether angles and thrust, and also 
minimizes the thrust-sum in the commonest cases of interest (at hover and for iden- 
tical helicopters with equal load sharing). 

We remark that the thrust-sum varies with three independent configuration 
variables (e^, 8 ^., 6 ^ 3 ) » and that it can be minimized with respect to all three 
variables. However, it is demonstrated below that the minimum given with respect 
to alone is at the absolute lower bound given by the apparent load-sum in the 

commonest cases of interest and otherwise is within 1 /t of this lower bound for a 
large set of empirical examples. Thus, and Bj. can be selected as desired with 
no significant effect on thrust requirements. 

In the analysis below it is convenient to minimize with respect to B, but this 
is equivalent to using either tether angle or F34 as the independent variable. 


A Suboptimum Configuration 


A near-optimum choice of B with simple results is obtained if we select 
B = 0. In this case, equations (7.2) and (7.7) give the tether and thrust forces 
as F13 ” . F24°, Tl^®, and T2®. These are repeated here for convenience. 

if: B = 0, then: F34 = -A 


FI 3 u13 = FI 30 = 0.5 mjj£3 + y m^^ft 

F24 u24 = F240 = 0.5 mj^fi* + (1 - 

XI = Xl® = -0.5 mj^f3 - V m^^n - m^n 

X2 = X2° = -0.5 m^^f4 - (1 - u)mj^fi - m^^fg 

where 

y = 0.5 cos(6 - 6j^)/cos 6 cos 

A = 0.5 cos(6 - t^)cos(S + e^)/cos cos 6 sin 6 


(8.3) 
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The results for any other choice of B are related to these by adding the term 
B along the spreader bar as given previously in equations (7.2) and (7.7). 

The results for this configuration are denoted by the superscript ( )® and are 
readily characterized: F34 depends only on the apparent suspended load but is 

independent of jQ, f2, f2» apparent loads appear linearly in the tether 

and thrust forces; the suspended load is distributed to the tethers and the helicop- 
ters in the ratio ii/(1 - y), which depends only on and is identical to the 
bridle-cable force ratio (fig. 5.1); and the spreader-bar load is distributed very 
nearly equally to the two helicopters. The tether tensions and directions are 
independent of the helicopter apparent loads and are generally aligned nearly paral- 
lel to ft (except for the unloaded tether when is an extreme and y or 

(1 - y) is zero), and the thrust for each helicopter is independent of the other 
helicopter. 

As is shown below, this simple configuration also minimizes IT in two cases 
of common interest. More generally, its thrust sum exceeds the minimum. Recalling 
the lower bound given by the apparent-load-sum (fig. 7.1), then the minimum thrust- 
sum, denoted IT*, is bracketed by 


IL < IT* < iLo 


(8.4) 


where 


ZT° = T1° + T2° 

IL = nij^fa + mj^fb + m^fM + m ^f2 
IL = 1^1 


Minimum Thrust-Sum Configuration 


The thrust sum is given in figure 7.1 as a sum of positive functions of B, 
each with a single minimum: 

IT = [(T1° + B)^ + + [(T20 - B) + 

where T19, T29, a., and are defined in figure 7.1. The value of B that 
minimizes IT is (algebra omitted): 


Then, 


B* = (a^T20 - 02T19)/(a^ + 02 ) 


(8.5) 


IT* = [(T1° + T20)^ + (a^ + 
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This result is derived assuming that fM, f2 do not vary with tether direction or 
rotor thrust. This assumption is not strictly satisfied, since the helicopter 
aerodynamic forces vary with thrust, but these variations are negligible. 

Some functional properties of the minimum ET configuration can also be 
given. The minimum point for each thrust is 


min(TI) = a. at B = -T19 
F34 ^ ^ 


min(T2) = a_ at B = T29 
F34 ^ ^ 

and then 


eT* > 0^+02 and B* is between -T1^ and T2^ 

These properties indicate that B* lies between the minimum points for the two 
thrust functions, and, in general, neither T1 nor T2 is minimized at B*. 

Expressions for the corresponding values of thrust, tether directions, and 
other variables that depend on B are listed in figure 8.3. For the minimum ET 
configuration, B* and the other variables depend on all the apparent loads of the 
system. If B* = 0, then the optimum configuration is given by equations 8.3. This 
occurs if Tl^* and T2® have identical tilt angles relative to the spreader bar 
(when T19/a^ = T 29 / 02 ); for example, when (1) ej^ = 0 and ra^n = m 2 f 2 (identical 
helicopters, equal loading) or (2) the apparent loads fl, n_, f2, f3, f4 are all 
parallel (e.g., at hover). 


The general results in figure 8.3 indicate that the optimum thrust angles 
e!|,e* are equal (equal tilt from the spreader bar); however, the thrust vectors are 
not necessarily parallel, since their angles out-of-the-plane of the triangle, 

X^,X 2 , can differ. In addition. 


I tan e!|| = 


lEL. 


+ a. 


|ELJ 

„,2J/2 

(eLj + ELj^) 


= I tan 


and then 


|e*l < UlI 

This inequality follows from the definitions of 0 ^ 02 * and the identity, 

T1 ° + T2<J = -EL; the result states that (e![| is less than the tilt of ^ from 
the spreader bar, e^. The thrust ratio p| depends on and the apparent loads, 
and can be controlled by choosing appropriately. 

These results are difficult to characterize with respect to several properties 
of interest, such as (1) the optimum tether directions, and (2) the effects of load, 
spreader-bar and helicopter aerodynamics, and apparent loads on the solutions. 
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DEFINITIONS:®'** 

^ = m£f2 + mbfb + mif1+m2!2 


SPREADER-BAR COMPRESSION: 

B* = T2p - «2 T1p)/(ai + ^2) 


B* = 0.5 mg ffi sin 5 


a-j + «2 


- tan 6 tan eg 1 - 0.5 fbj 


a.| - «2 a.j m2 f2j - 0:2 m>j fij 


“2-“l 




+ a<i 


F34* = -A + B* 


THRUST FORCES®: 

TlJ = T1? + B*(1,0, 0)^ = 

T2j =T2^ -B*(1,0,0)''’ = 


“1 ^*-1 
a-| + «2 


«2^*-i 

^2 



tan e!| = tan e| = -SLj/(a.j + 02) 

T1* = ».|/co$ e^ 

T2* = 0 : 2 / 00 $ e^ 

2 i2i 

2T*= [SL> + («! +a2rl 
Pj* = 0 £.j/o:2 


TETHER FORCES:® 

F13* u13* = F13^ - B*(1, 0, 0)'*' = 

F24* u24^ = F24? + B*(1, 0, O)'*' = 

tan e^3 = (a>j tan e 3 | + m-| flj)/oi3 
tan C24 = («2 ^2 '*’2 ^2j)/a4 

F13* = 0:3/00$ eJ|3 
F24* = a4/co$ e^ 


\“l+“2 


\ “1 + “2 


-m., fij. F13j*, F13g 
- m 2 f2j. F24f, F24g 


® M. A. F13°, F24°, H® 11°, ^2, 03, 04 ARE DEFINED IN FIG. 7.1. 

ASSUMPTIONS LISTED IN FIG. 7.1 ARE MADE HERE. 

® X^. X3, X.|3. X24 ARE INDEPENDENT OF €13 AND ARE GIVEN IN FIG. 7.1. 

Figure 8.3.- Minimum thrust-sum configuration: tether and thrust forces. 
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However, two special cases for which the minimum ZT configuration is given by 
equation (8.3) are more easily interpreted and are considered next. 


Results for Identical Helicopters and Equal Load Sharing 

It is expected that the most common dual-lift case will be the one in which the 
helicopters are identical and there is equal load sharing. Results can be derived 
from the equations given in figure 8.3 using the following assumptions: 

“l = ”2 = \ 
fi = f2 = £Ti 

f3 = fi = fb 

The equal helicopter masses are denoted m^j. Equal apparent helicopter loads, 
denoted fh, can be assumed for identical helicopters with equal loading, but it 
suffices that their difference be negligible ( |fM - f2| << 0.1 g). The differ- 
ences |f3 - fb| = |f4 - ^1 = 1^1 are zero in static equilibrium and generally 
negligible over the maneuvering domain of interest. Last, the spreader bar is set 
level relative to the apparent suspended load = 0), and this results in equal 
helicopter loading (p = 0.5, T1* = T2*). 

These assumptions result in considerable simplification of the equations given 
in figure 8.3 because = a^, and a_ = a^,. In addition, B* = 0, so that the 
optimum configuration is given by equation (8.3) for this case. Results for the 
tether and thrust force variables are given in figure 8.4. Results are also given 
in terms of the dimensionless apparent-load difference vectors, in order to 

indicate the effects of aerodynamic forces, where 

Itb = m|^(fb - fll)/mj^ft = mj^(fab - fat + ^ - ab)/m^ft 

th = - ft)/mj^ft = m^ ( fah - fat + at - ah)/m^ft 

These differences are due to differences in aerodynamic specific forces and acceler- 
ations, and, in view of the discussion in section 4, are negligible at hover and low 
speeds, and are dominated by aerodynamics at higher speeds (principally as drag 
differences). As shown below, the tether and thrust directions depart from the 
direction of ft in proportion to these differences. Note that for m^^ << m^^, then 

I tb| << 1 in all circumstances within the assumptions of figure 7.1. The following 
properties are seen in figure 8.4(a): 

1. The optimum tether directions are parallel; each tether carries half the 
apparent load attributable to the suspended load and spreader bar and is parallel to 


47 


ASSUMPTIONS: 
eg = 0 
nil = ni2 = 

f1=f2 = fh OR ln-f2l«:0.1g 
f3 = f4 = fb OR lAfl«0.1g 


DEFINITIONS: 

-b “ ~ 

Cj, = - fe)/mgfe 


RESULTS: 

IXpl, le*l<7r/2 ne {1.2. 13, 24} 
F34* = -O.BmgfCCtnS 


tan Xi-j = tan Xo.* = 


nibfbj 


Cb: 


24 1 m|j/mg + Cb|^ 


tan 6^3 = tan €24 ” 




-Kb: 


y/ (m|jfbj)2 + (mgfC + ni(jfb|^)2 -y/ Cb? + (1 + nrijj/mg + Kb|^)^ 

F13* = = 0.5(mgfC + nijjfb) 


tan Xi = tan Xo = 


nibfbj 2m^,fhj 


Cbj + 2Chj 


I 3 mgfC + + 2m|^fh|^ 1 + nij^/mg + 2mj^/mg + Kb|^ + 2Ch|^ 

-nifatbi - 2mbfhj 


tan ei = tan e% = 


\J (nibfbj + 2nibfhj)^ + (nigfC + nibfbj^ + 

-Cbj-2Chj 

x/(Cbj + 2Chj)2 + (1 + mb/mg + 2nib/mg + Cb,^ + 2Ch,^)2 
T1 * = 12* = -0.5 mgfC - 0.5 tiibfb - ^htil 
2T* = ImgfC + mbt^+ 2nibfhl 
P^-=1 


(a) EQUATION SUMMARY 


Figure 8.4.- Minimum thrust-sum configuration: identical helicopters, equal load 

sharing. 
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LOAD DRAG, d, g 


(b) EFFECT OF LOAD DRAG ON THRUST DIRECTION ANGLES: 
(eg, 6, m|y mj,) = (0, 60°, O.lmg, l.Bmg) 

fab = fah = ao = 0 
fag = -d 


Figure 8.4.- Concluded 
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this load vector-sum. This tether direction is independent of the helicopter 
apparent loads. The tether angles, ^*3* given in terms of the 

spreader-bar load-difference vector, Since |^| << 1, these angles are small 

(approximately tb., ab^, respectively) and the tethers are close to k^. and ua 
(nearly in the triangle plane and nearly perpendicular to the spreader bar). These 
angles increase with the spreader-bar mass ratio mj^/mj^. 

2. The required helicopter thrust vectors are parallel; each equals half the 
total apparent load and is directed opposite this vector-sum. The thrust-sum equals 
the lower bound given by the magnitude of the total apparent loads so that 
ZL = = Z^ for this case. The thrust direction angles relative to triangle axes, 

X^, e!|, are also given in terms of the difference vectors, tb^, J,h^. These angles 
are small, and the thrust direction is nearly aligned with provided 

1^1 << 1. However, moderate differences in aerodynamics, |fah - f^|/f4, can occur 
and can result in larger angles between k^. and u1*. u2*. This is illustrated in 
figure 8.4(b) by the effect of load drag on these angles for various values of Sf 


Results for Parallel Apparent Loads 

For the parallel apparent-loads case the apparent loads of all bodies in the 
system are parallel; that is. 

Parallel apparent loads: £1 = f1 ut, f2 = f2 ua, f3 = f3 ui, f4 = f4 u£ (8.6a) 

Windless hover is a case of general interest which satisfies equation (8.6a) and for 
which the specific apparent loads are also equal in magnitude: 

Windless hover: fat - fab = fal = fa2 = 0 

cjv = 0 

(8.6b) 

a& = ab = a1 = a2 = 0 

and then 

fa = fl=f2 = f3 = f4 = gkjj 

For parallel apparent loads, B* = 0 and the optimum configuration is again 
given by equation (8.3). Results for all variables that depend on B are listed in 
figure 8.5, including values at the extremes of and at = 0. In addition, 
graphical results are given figure 8.6. These show the following properties of the 
minimum Zq. configuration for this case: 

1. The tethers and thrust vectors are aligned parallel to the apparent sus- 
pended load vector, fa, at all This is always in the plane of the triangle and 

at an angle from the perpendicular to the spreader bar, and the system is 

coplanar . 
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ASSUMPTIONS: 


f1 = f 1 u«, « - f2 u«, f3 = f3 u«, f4 = f4 uC 
DIRECTIONS: 

u13 * = u24* - u1* = u2* = u« 

«1*3=«24“ «1“"2 = «)2 
^13 “ ^24 = ' ^2 “ ® 


FORCES: 


SYMBOL 

GENERAL FORM 

VALUES AT 

eg = -7 t/2 + 6 

ej = 0 

eg « ff/2 - 6 

F34* 

COS(6 - € q ) C0S(6 + €n) 

0.5 mp ffi i ^ 

^ cos 6 cos eg 

0 

-0.5 mg ffi ctn 6 

0 

M 

0.5 cos(6 - eg)/cos 8 cos eg 

0 

0.5 

1 

F13* 

ju mg + 0.5 f3 

0.5 mjj f3 

0.5(mg ffi + mjj f3) 

mg ffi + 0.5 mjj f3 

F24* 

(1 - fi) mg ffi + 0.5 mjj f4 

mg ffi + 0.5 m|j f4 

0.5(mg ffi + mjj f4) 

0.5 m^ f4 

T1* 

jLt mg + 0.5 m|j f3 + m-| f1 

0.5 m|j f3 + m-| f1 

0.5(mg ffi + m^j f3) + m^ f1 

mg ffi + 0.5 mjj f3 + m.^ f1 

T2* 

1 {1 - /i) mg + 0.5 mjj f4 + m 2 f2 

mg ffi + 0.5 mjj f4 + m 2 f2 

0.5(mg ffi + mj, f4) + m-| f2 

0.5 m|j f4 + m 2 f2 

ST 

mg + m^j fb + m^ f 1 + m 2 f2 





Figure 8.5.- Minimum thrust-sum configuration: parallel apparent loads. 

2. The spreader bar is in compression at all le^^l < fr/2 - 6. The compression 
is zero at the extremes of and maximum at = 0. This maximum is propor- 
tional to ctn 6 and decreases with 6. 

3. Each tether supports half the spreader-bar apparent load plus a fraction of 

the suspended load determined by At the extremes of ej^, one tether supports 

the entire load, and at = 0 the tethers support the load equally. 

4. Each helicopter supports its own apparent load plus that portion of the 
spreader bar and suspended apparent loads carried by its tether cable. 

5. The thrust-sum equals the sum of apparent loads, m fa + m. fb + m.fl + m,f2, 

at all e^. a b i 

Thus, for parallel apparent loads, the optimum geometry is simple; the entire 
system is coplanar, with tether- and thrust-force vectors parallel to the apparent 
loads, and with thrust requirements at the lower bound, ZL. 
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(a) SPREADER BAR FORCE 



(b) TETHER FORCES 


— F13* + F24* 

^ + mjj)fC 

5, deg 
■^45 

^ F13* 

“ (mj + mb/2)f« 


6, deg 
/45 
. .50 


(mg + m|j/2)fS! 


+ T2* 

^ + m2 + trig + m|j)fC 

deg 

T-i* 

!^x:r^50 

(m^ + mg + mjj/2)fg 

60 

45 T2* 

55 (ff'? + mg + mL/2)fC 

.60 


eg, deg 

(c) HELICOPTER THRUST 


Figure 8.6.- Minimum thrust-sum configuration: results for equal specific apparent 

loads (f1 = f2 = f3 = = ft), (mjj,m^,m2) = (0. 1 , 1 .5, 1 .5)mj^, 6 = 60°. 



General Results for Minimum Thrust-Sum Configuration 

Parallel loads are obtained at hover and result in the simple geometry noted 
above. It can be shown that if the angles between these vectors are small, or if 
the apparent load is small compared to mj^fJ., then the properties given above for 
parallel apparent loads approximately characterize the minimum thrust-sum 
configuration. 

The spreader-bar apparent load is always small (mj^fb << mj^ft) in the present 
context so that significant departures from the results in figure 8.5 require sig- 
nificant angles between the apparent loads of the helicopters and cargo. In most 
circumstances, these angles are small in equilibrium flight. Referring to the 
sketch, the angle between f, ft satisfies 

If - fil I fa - fM + M - a| 

< = ri 

where f, fa, and a are associated with either helicopter. The differences in 
specific apparent loads are due to differences in aerodynamic forces and c.g. accel- 
erations. The discussion in section 4 indicates that acceleration differences are 
of the order of 0.1 g or less on all segments of the trajectory for the maneuver 
domain and system sizes of interest here, and the aerodynamic forces are dominated 
by drag and tend to be parallel so that differences, |fa - fan, | . of the order of 
0.1 g require sufficiently large airspeeds, and are unlikely to exceed 0.5 g at any 
speed. Thus, the properties noted in figure 8.5 provide a good approximate charac- 
terization of the minimum ZT configuration to moderate airspeeds, of the order of 
50 knots. 

Computations were made to indicate the variation in configuration geometry with 
apparent-load differences. The algorithm used to compute the configuration is given 
in figure 8.7 and comprises a solution of the restricted trim problem stated in 
section 3. The input quantities consist of system parameters, the formation angle, 
and spreader-bar tilt, and the specific aerodynamic forces and reference trajectory 
acceleration. Simplifying approximations are made in computing the apparent loads; 
that is, differences in c.g. accelerations of the constituent masses and the 
spreader-bar rotation effect ^ are neglected. Accurate calculation of these 
quantities and the input aerodynamic forces requires solution of the general trim 
problem outlined in section 3, but reasonable approximations can be made for the 
present illustrative calculations based on section 4. The remaining equations are 
taken from equations (6.3) and (6.4) and figures 7.1 and 8.3; these are valid for 
any equilibrium configuration except that the parameter B is computed to mini- 
mize zT. 

Cargo drag is the principal aerodynamic force in equilibrium, and results were 
obtained for its effects on the configuration; that is, for the case in which 
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INPUTS:® 

Ftig, tTijj, mi, m 2 , 5, eg, Pf 

fa£^, fab|i, 

SPECIFIC APPARENT LOADS, LEVEL-HEADING AXES:** 

“ 9H " ®°H 

f3H = fb„ 
f4H = fbH 

flH=gH + falH-®°H 
^2 h “ 9H " ®°H 

TRANSFORMATION OF TRIANGLE AXES: 

fe=l«Hl 

uIh = 

= sin~^ 

2 = tan~^ I(u®x 

= 2 + sin~^ (sin eg cos 2/u£2) 1^1 - 2 | < jr/2 

■•■t.H = El (0t> E2 E3 

SPECIFIC APPARENT LOADS, TRIANGLE AXES: 
f£| = f£ {-sin eg, 0 , cos eg) 

H 

AUXILLIARY QUANTITIES: 

ti = 0.5 cos (6 - eg)/cos 8 cos eg 
F13^ = 0.5 m|j f3i + /img f£i 
F24° = 0.5 mb fA^ + (1 - m) mg f£^ 

T1^ = -F13^-mi fit 
T2? = -F24^ - m 2 

Figure 8.7.- Restricted trim algorithm: minimum z.j. configuration. 
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2 2 

ai = (Tlj» +Tlg ) 

2 2’-^ 

«2 = (T2J» + T2g ) 

BRIDLE-CABLE TENSIONS: 

F35 = 0.5 mg fC cos(5 - 6g)/co$ 5 sin 6 
F45 “ 0.5 mg fH cos(S + eg)/cos 6 sin 6 

SPREADER-BAR COMPRESSION, MINIMUM ST CONFIGURATION: 

B* = (a, T2f - «2 Tif )/(a, + «2) 

B = B* 

F34 - B - 0.5 mg f£ cos(5 - eg) cos(5 6g)/cos eg cos S sin S 

THRUST AND TETHER FORCES:® 

T1t = T1° + B(1,0,0) 

T2t = T2g-B(1,0, 0) 

F13t=F13g-B(1,0, 0) 

F24^ = F24° -I- B{1, 0, 0) 

® ASSUMPTIONS LISTED IN SECTION 2 ARE MADE HERE. 

^ DIFFERENCES IN C.G. ACCELERATIONS AND Af ARE NEGLECTED. 

® ^ ARE TETHER FORCES AS APPLIED TO SPREADER BAR. 

Figure 8.7.- Concluded. 

fan = -d ij^ 

fab = fa1 = fa2 = 0 

where drag is represented as a force along the ground-track direction for these 
computations (this assumes zero wind and y for simplicity). Computations were 
made for 

d e [0, 0.6] g 

le^l < ir/2 - 6 

e [00, 900] 
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and for both static equilibrium and turns (a = 0 and a = 0.3 g for all 
bodies), and for both identical and nonidentical helicopters (m< = m 2 = 1.5 njj^ 
and m^ = m^^, m 2 = 2mj^). Some results are presented in figure 0.8 for the case of 
identical helicopters, static equilibrium, and s^. = 45°. 

Results for tether- and thrust-direction angles (fig. 8.8(a)) show that 

increases with cargo drag at all but is small (under 5°) except when the 

tether is unloaded (ej^ = -30®); that is, tether No. 1 is always nearly in the tri- 
angle plane except in the extreme case in which it is unloaded. The thrust angle 
|xj increases with cargo drag to moderate angles (15® - 20®) at all so that 

the thrust vector is more noticeably inclined from the triangle plane as drag 
increases. The tether and thrust tilt angles, ^tre shown versus e^. 

For zero drag, these directions are parallel to ^ (e-l 2 * = = Gj^) as expected; 

for nonzero drag these angles depart from this simple relationship increasingly with 
drag. This departure becomes large near the extreme case in which the tether is 
unloaded. These conclusions apply to both tethers and thrust vectors in all cases 
calculated. 

Additional results in figure 8.8 show that the angles between the two tethers 
and the two thrust vectors, z.( u13 *. u24 *) . z.(T1*, T2*) , are zero for (1) = 0 

independent of drag and (2) d = 0 independent of Gj^. That is, these are parallel 
vector pairs for the cases of ( 1 ) identical helicopters with equal load sharing and 
(2) parallel apparent loads, as previously noted in figures 8.4 and 8.5. In other 
cases, these are not parallel, and moderate mutual tether and thrust angles must be 
maintained by the helicopters as a function of drag, helicopter masses, 0^ and g^^. 

The angle between Tl^* and apparent-load-sum was computed (not shown). The 
thrust vectors are necessarily along whenever ZT* is at the lower bound, 

zL. This angle was found to be zero for = 0 or d = 0 and otherwise increases 
with drag but is under 10® in all cases calculated. In addition, the required 
thrust-sum and its lower and upper bounds (eq. (8.4)) were also computed. These 
were within of each other in all cases computed; that is, both ZT* and ZT® were 

within of the absolute lower bound on thrust requirements given by ZL. 

In general, ZT* varies with the independent configuration variables not 
included in the optimization, Sj and Gj^, but the present empirical results and 
special-case analytical results above demonstrate that the variation is weak or 
nonexistent. The possibility of variations in the aerodynamics and, hence, in the 
apparent loads, fi, fb, f^, and f2, with 6^., Gj^ is not considered in the preceding 
statement since the apparent loads enter the problem as parameters in this report. 
However, it can be assumed that these variations are negligible for the load and 
helicopters and that the spreader-bar aerodynamic force may vary strongly with 6^ 
but always satisfies I Fab I /W^ << 1; therefore, its variations are negligible in the 
sums ZL,ZT*. Thus, zT* is very nearly invariant with 6^ and Gj^, and no signifi- 
cant further reductions are available by minimizing ZT* with respect to these 
variables. 
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9. LOAD DISTRIBUTION 


In this section, load distribution refers to the apportioning of the thrust 
requirement IT* to the two helicopters. This can be controlled within a range of 
values by selecting the spreader-bar tilt, appropriately within its admissible 
range for no cable collapse. However, each helicopter has a thrust limit, Tl^^jj, 
T2max» which further restricts the range of load distributions that can be flown to 
those for which all thrust margins are positive; that is, can be flown provided 
that 


where 


M,U,) 

Kj(ej) 


T1. 


max 


T2 


"max 


IT 


max 


- T1*(e,) 

- T2»(e^) 

- lT*(e^) 


> 0 
> 0 
> 0 


IT 


max 


A 


T1 + T2 
max max 


(9.1) 


To examine the corresponding limits on load sharing, consider the space of 
thrust combinations (see fig. 9.1) and let p be the locus of thrust combinations 
for the minimum IT configuration generated as varies over its admissible 
range : 


p = [(T1*,T2*): -if/2 + « < e- < if/2 - 6] 

Xr 


(9.2) 



Figure 9.1.- Available and required thrust combinations. 
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Assume two properties of p: 

A1. T1* increases with 
A2. T2* decreases with ej^ 

The endpoints of this curve-segment, p^, P2» correspond to the extremes of 


e 


E 


-ir/2 + 6 
ir/2 - 6 


Next, consider the helicopter thrust limits 
thrust combinations given by (see fig. 9.1): 


T1 


max’ 


T2, 


max’ 


and the region of 


9 - - (0. X [0, T2^^] 


Every combination in this region satisfies the requirements for positive thrust 
margins (eq. (9.1)), and every point outside does not. Therefore, only that subseg- 
ment of p inside this region, if any, can be used in selecting Ej^. Assume that 
the use of dual lift is restricted to suspended loads that cannot be carried by a 
single helicopter. Then, T1 „,^ y, T2 „,^y must be less than the maximum values of 
T1*, T2* at the endpoints of p since these correspond to unloading one of the 
helicopters. In addition, assume that dual lift is possible; that is, that p 
intersect the required region. For this, it suffices that ZT* be less than 
zTmax This requires that p be below the boundary T1 + T2 = 

in figure 9.1. This is not an unduly restrictive assumption since IT* is very 
nearly invariant with Ej^ for all the special and empirical cases discussed in 
section 8. Thus, for situations in which dual lift is both necessary and possible, 
assume 


A3. < T1»(ej^), T2|^^jj < T2*(£j^) (dual lift is necessary) 

A4. zT*(ej^) < for Ej^^ ^ (dual lift is possible) 

The assumptions, A1 - A4, suffice for p to have some nontrivial subsegment within 
the required region: 


p' = [(T1»,T2»): E < E < E ] 

where the endpoints of this segment, p^, Pij (see fig. 9.1) correspond to tilt 
angles, e , e , defined by 

»3 


(9.3a) 


T2*(e, ) 
3 


T2_ 


max 


(9.3b) 


T1»(e. ) 


T1 

max 
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Then the available range of the thrust ratio, T1*/T2*, is 


Ue. ) 


IT» 


T1. 


T2 


max 


- 1 < p| < ) = 


max 


ZT* - T2 


max 


(9.3c) 


A suitable choice of thrust ratio is given by the thrust limits: 

= T1 /T2 
Tx max max 


(9.M) 


This occurs at corresponding to the intersection of the locus, p, with the 

line T1/T2 = figure 9.1, and it is always a possible choice (it is 

always in the subsegment, p' ) for the present assumptions. Note that this choice is 
also mandatory when the load equals the system's maximum load-carrying capacity 
(when IT* = since the only such point in the region of non-negative thrust 

margins in figure 9.1 occurs at (Tljjj^^j^, Thus, only when the suspended load 

is significantly below the system's maximum payload is there any significant range 
of choice of p| other than p.j.jj. The thrust margins for this choice are in the 
same ratio; if 


if: p» = p^j^ , then: = p^^ 

The conditions A1 - A4 suffice for a solution to exist that has positive 
thrust margins. The results for parallel apparent loads (figs. 8.5 and 8.6) sat- 
isfy A1, A2, assuming that the apparent loads are constant or that they vary 
negligibly with e^. In addition, ZT* = ZL and is invariant with for this 

case; the locus p is the straight line T1 + T2 s zL; and the condition A4 
requires simply ZL < ^ solution to exist. 

The conditions A1 and A2 and the insensitivity of ZT* to are difficult to 
verify analytically for arbitrary apparent loads. However, these conditions are 
satisfied by all the empirical cases discussed in section 8, and at flight condi- 
tions in which the apparent loads are nearly parallel, which encompasses low-speed 
flight and is expected to encompass most flight conditions and helicopter-load 
combinations of interest. 


Further, it can be shown that if the configuration is selected to obtain 
B = 0, then T1® and T2® are, respectively, increasing and decreasing with ej^. 
Thus, conditions analogous to A3 and A4 suffice for solutions with positive thrust 
margins for this configuration also. 

Finally, assuming that selected to minimize ZT, then the value of 

corresponding to any choice of p| can be calculated from the general equations 
in figure 8.3: 
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p? = — 




T 1/2 


O.Snjj^fft cos(« 


e^)/co3 «]‘ 


^(O.Smj^fMj + m 2 f 2 j)^ + [0.5mj^f4j^ + ®2^^k * cos(6 + £^)/cos 


( 9 . 5 ) 


Empirical inversion appears necessary in general (note that the triangle axes compo- 
nents of £2, etc., vary with e^^), but a closed-form solution can be given for the 
case of parallel apparent loads from the results in figure 8.5: 


tan e. 
i 


ctn 6 


p* - 1 m2f2 - m^fl + ni^^Af 

+ 1 m^^ft m^ft 
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10. GENERAL EQUILIBIUM CONFIGURATIONS 


The general equilibrium results for arbitrary tether angles given in figure 7.1 
are reconsidered in this section to determine (1) the variation in system geometry 
and thrust requirements with the tether tilt angle ^2) the sensitivity of 

the geometry and thrust requirements to variations in the tilt angle from its opti- 
mum (minimum ZT) value. The previous results can be rearranged in terms of values 
at the optimum configuration and the effect of variations of the No. 1 tether tilt 
angle. Results are listed in figure 10.1 along with derivatives to indicate local 
sensitivity of the equilibrium configuration. The results are also illustrated 
graphically in figure 10.2 for an example case with unequal helicopters and drag 
forces on all bodies. 

The general variations in tether force and thrust vectors with tether angle are 
readily characterized from the analytical results (fig. 10.1); only their components 
along the spreader bar vary with tether angle, and these variations are identical in 
magnitude to the variation in spreader-bar force for all four forces with opposite 
signs for the two tether forces and the two thrust vectors. In addition, the fol- 
lowing properties of the variations can be seen. 

First, the spreader-bar force increases linearly with the tilt-angle variation, 
expressed as tan - tan e^_. Its proprotionality factor is an increasing 

function of that is, sensitivity to variation of the No. 1 tether tilt 

increases as the tether carries more of the cargo in the optimum configuration. 

At = -it/2 + 6, the No. 1 tether is unloaded and << m^^ft, whereas at 

= if/2 - 6 the No. 1 tether carries the entire load and sensitivity increases an 
order of magnitude to o, “ variation in sensitivity with is 

evident in figure 10.2(a7 and also occurs for the other variables (e.g., 
figs. 10.2(b) and 10.2(d)). When expressed in terms of tilt-angle variations of the 
No. 2 tether, tan - tan (fig. 10.1), F3^ is found to be linear in this 

variation with proportionality factor, which decreases with ej^; that is, sensi- 
tivity to €24 increases with increased loading of the No. 2 tether. 

The results (figs. 10.1 and 10.2(a)) also indicate that the spreader bar can be 
placed in tension for sufficiently large tilt-angle variations and that F34 
becomes arbitrarily large at any as |ei3l approaches it/2. However, the 
realizable range of F34, tether angles, and system forces in equilibrium is limited 
by the helicopters' maximum thrust capabilities. 

Second, the No. 2 tether angle, tan 624, varies linearly with tan 
proportionality factor, -00/04, is negative [that is, ®24^®13^ 

opposite directions from the optimum configuration (e*o» ®24^®13^^^ increases 
with At = -ii/2 + 6 the No. 1 tether is nearly unloaaed, o^/o^ << 1, and 

even large variations of require very little change in 624 (fig. 10.2 (b)). 

However, at - n/2 - 6, the sensitivity is very large, 0^/04 >> 1, and the 
required angle |e24l driven to values near v/2 by even small variations, 

e - E* 

®13 13* 
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VARIABLE 

GENERAL FORMULA 

DERIVATIVES^ ^ 

\d tan £ 13 / 

F13u13 

F13* u13* - (F34 - F34*) 

-“3 it 

F24 u24 

F24* u24* + (F34 - F34*) 

“3 it 

T1 

I 1 * + (F34-F34*)1i 

“3 it 

12 

T2*-(F34-F34*)ii 

-«3it 

F34 

F34* + tt 3 (tan £13 - tan £ 13 ) 
F34* + tt 4 (tan £34 - tan £ 34 ) 

“3 

tan £24 

ao 

tan £34 - — (tan £13 - tan £ 13 ) 

“4 

F13 

F13* cos £i 3 /cos £i 3 

«3 sin ^13 

F24 

F24* cos £24/cos £34 

-«3 S'" ^24 

2F 

F13+F24 

a 3 (sin £ 13 * - sin 634 *) 

tan £1 

^3 

tan £1 + — (tan £13 - tan £ 13 ) 

“1 

tan £2 

Oio 

tan ~ ^13 " ^ 13 ) 

'“2 

T1 

T1* cos £i/cos £1 

tt 3 sin £1 

T2 

T2* cos £i/cos £3 

-Q !3 sin £1 

2T 

T 1 +T 2 

0 


DEFINITIONS: 

«1 = 0-5{Kf3j + mi flj)^ + [m|j f3|^ + mi + mg f8 cos(6 - eg)/cos 5] ^ | ^ 
a 2 - 0.5 1 (mjj f4j + m 2 f2j)^ + [mjj f4|^ + m 2 f2|^ + mg f£ cos(6 + 6g)/cos 6 ] 

«3 = 0.5 1 (mjj f3j)^ + {m|j f3|^ + mg f£ cos(6 - 6g)/cos 

= 0.5 I (m|j f4j)^ + [mjj f4|^ + mg f£ cos (6 + eg)/cos 5] ^ ^ 

Figure 10.1.- Equilibrium configuration with tethers misaligned from the direction 

for minimum ET. 
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F24 F13 

tC + 0.5mhfbl ImnfC + O.Bmufbl 











(f) HELICOPTER THRUST SENSITIVITY; TETHER 
MISALIGNMENT FOR 1% PENALTY 

Figure 10.2.- Concluded. 

If given in terms of the misalignment of the No. 2 tether, 


tan e 


13 


tan “ 


(tan 


■’2H 


- tan 


it is seen that the reciprocal proportionality and the reverse sensitivity trend 
with occurs. Thus, near the extremes of load sharing, the lightly loaded 
tether is required to undergo large variations in tether angle between ±90° to 
maintain equilibrium in the presence of even small changes in the direction of the 
heavily loaded tether. 
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Third, the general variations of tether tensions are more easily anticipated 
from their relations with F34: 

? ? 1 /p 

F13 = [(F34 + a - F13^) + “ 3 I 
F24 = [(F34 + a + F249)^ + 

These are positive functions of F34 with a single minimum point and with arbitrar- 
ily large values as |F34| -► ® or, equivalently, as behavior is 

seen in figure 10.2(c). The line of optimum configurations is included in the 
figure. It coincides with the minimum points of F13 and F24 when their deriva- 
tives are zero (fig. 10 . 1 ); this occurs at = ^24 " ^ at hover with 

= 0) and is very nearly obtained when these angles are both small (identical 
helicopters with = 0). More generally, these derivatives have opposite signs 
whenever ^* 3 >^24 same signs, and (locally) one tension increases and the 

other decreases with " ^* 3 * 

The derivatives of the tether-tension-sum (fig. 10.1) are zero and the tether- 
tension-sum is minimized at the optimum configuration simultaneously with ET 
whenever occurs at hover and for identical helicopters with 

H = 0. 

Fourth, the thrust angles, tan tan t2* vary linearly with the tan 
in mutually opposite directions (fig. 10.1). Their sensitivities to tilt variations 
(fig. 10 . 2 (d)) increase with Ej^ and are both smaller at nearly all Ej^ than the 
sensitivity of E 24 (compare figs. 10 . 2 (d) and 10 . 2 (b)) because of the influence of 
the helicopter apparent loads in defining the required thrust-vector directions. 

Last, the helicopter thrust magnitudes are positive functions with a single 
minimum and arbitrarily large values at 1 ^ 13 ! = 

figure 10.2(e) for the example case. In general, the derivatives of T1 and T2 are 
equal and opposite at the optimum configuration (fig. 10 . 1 ) so that locally one 
thrust Increases and the other decreases. This occurs because the optimum eI^^ 
always falls between the values of e^^i which minimize T1 and T2 separately. 

These derivatives are zero, and T1 and T2 are minimized, as well as ET, 
whenever e!^ = 0 or, equivalently, when the spreader bar is perpendicular to 
EL (when ELj^ = 0); this occurs at Ej^ = 4° for the example case, and usually 
occurs at small Ej^. These derivatives are small (relative to pei' radian) 

when e!| is small (when |eL^| << (a^ + 02 )). 

The thrust-sum derivative is zero at the optimum configuration by construction, 
and its sensitivity to tilt variations is small in the example (fig. 10 . 2 (e)). 
Sensitivity is examined further in figure 10 . 2 (f), which shows the boundary of 
misalignments corresponding to a thrust-sum penalty of IJt. This encloses a region 
that narrows with e^^ (from ± 25 ° at = 0 to - 18 °, + 12 ° at e,^ = 30 °) and indi- 
cates a moderate increase in sensitivity with increased loading of the No. 1 
tether. Similar calculations were made for each helicopter thrust, and the boundary 
for which both variations are under of the optimum thrust is included in 
figure 10.2(f). This boundary is inside the corresponding boundary for ET at 
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all 6j^; it meets that boundary in the vicinity of = 4® where the derivatives 
of T1 and T2 are zero, and it narrows considerably to -3°, +2® at = 30®. 

These results indicate a broad optimum with respect to the choice of tether 
angle at and near equal load sharing and that the thrust requirements (T1, T2, ZT) 
are not sensitive to tether-tilt variations, £12 “ ^*3 ®24 ~ ^24 this 

case. Thus, the configuration corresponding to B = 0 with simpler formulas than 
the optimum (eq. (8.3); fig. 8.3) can be used for system flight coordination with 
little penalty. More generally, it is found that the sensitivity of system geometry 
and forces to misalignment of a tether from its optimum direction varies in propor- 
tion to the fraction of the apparent load that it carries, and that sensitivity to 
one or the other tether angles, e^2 ®24» increases with unequal load sharing. 

Near the extremes of e^, the equilibrium direction of the unloaded tether becomes 
very sensitive to variations in the direction of the loaded tether. This suggests 
difficulty in flying a system if the load distribution differs by as much as an 
order of magnitude. 
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1 1 . CONCLUSIONS 


The characteristics of a dual-lift helicopter system in equilibrium flight 
along any reference trajectory within its operational domain were examined. The 
system consists of two helicopters with a spreader-bar suspended below the helicop- 
ters on cables and the cargo suspended by cables from the ends of the spreader-bar. 

System characteristics (suspension geometry and forces, thrust-vector require- 
ments) are obtained by solving the force-balance equations at various points on the 
configuration, assuming that the c.g. accelerations and aerodynamic forces on the 
component bodies of the system are known. It can be shown that these quantities can 
be determined exactly, or to a good approximation, from the reference trajectory and 
from the force and moment descriptions of the bodies independent of solving the 
force-balance equations. Closed-form solutions of the force-balance equations can 
then be given with the c.g. accelerations and aerodynamic forces appearing paramet- 
rically. For each body, these quantities appear only as the combination of weight, 
inertia reaction, and aerodynamic force, which in this report is termed the apparent 
load (mf = m(g - a) + Fa). 

The detailed analytical results are summarized in the restricted trim algorithm 
(shown in fig. 8.7) and in the following conclusions: 

1. At any given flight condition, the force-balance equations have multiple 
solutions corresponding to three redundant configuration angles which can be 
selected by the pilot or autopilot. These can be taken as the spreader-bar heading 
relative to the ground track (formation angle, B^.); the spreader-bar tilt relative 
to the apparent suspended load vector (ej^), which determines load distribution; and 
the tilt of the No. 1 helicopter's tether relative to the spreader bar (£13) ♦ which 
determines the thrust-sum requirements. In equilibrium flight along a reference 
trajectory, the relative positions of the system's bodies are fixed or very nearly 
fixed (relative to axes attached to the load-spreader bar triangle) on any quasi- 
steady flight segment, with movement to new relative positions during transient 
maneuvering between these segments. This assumes the redundant configuration varia- 
bles are maintained invariant on these segments. 

2. Some limits on the parameters of the problem are required for the existence 

of solutions with appropriate properties. First, one bridle cable becomes slack if 
the other is aligned with the apparent suspended load, and both bridle cables are in 
tension only if the tilt |ej^| is less than tt/2 - 6, where 6 is the angle 
between the spreader bar and either bridle cable. Second, the triangle is right 
side up at any B^ and at any admissible if the angle between the apparent 

suspended load and the vertical is less than 6. This limits the combined c.g. 
acceleration and aerodynamic specific force of the suspended load for which the 
suspension is appropriate. Third, the tethers are in tension with the helicopters 
on the opposite side of the spreader bar from the cargo and their thrust is directed 
away from the spreader bar if the specific apparent load differences between the 
cargo and the other bodies are bounded; that is, if |fl_ - ft|, |f 2 - ft|, 

Ifb - ft| < ft sin 6 . 
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3. The inertial orientation of the triangle formed by the spreader bar and the 
suspended load varies along a reference trajectory and is defined by the condition 
that it always contain the apparent suspended load at the selected values of 

It is expected that the system will be flown with a fixed or slowly varying 
formation angle. The required spreader-bar pitch and load-swing variations with 
maneuvering and load aerodynamic force depend principally on 6^.; if the formation 
is longitudinal or transverse to the ground track, then turns or speed-change maneu- 
vers and lateral or longitudinal aerodynamic forces correspond to decoupled varia- 
tions in bar pitch or load swing from their attitude in unaccelerated low-speed 
flight. For intermediate formation angle, the required load swing and bar tilt are 
coupled . 

4. The required thrust-sum depends on the tether angle and is very nearly 
invariant with and Sf The tether angle can be selected for minimum thrust-sum; 
this maximizes the system thrust margin and also minimizes fuel rate at most flight 
conditions in the case of identical helicopters. A lower bound on the required 
thrust-sum is given by the magnitude of the sum of apparent loads; 

ZL = Inij^fa + m ^^fb + m.fj[ + m^f2 1 . This value is reached for special cases of inter- 
est (for identical helicopters with equal loading and at hover) and the correspond- 
ing thrust vectors are along -^. More generally, in empirical results for a 
comprehensive set of cases, the minimum required thrust-sum was found to be within 
1^ of this lower bound with small angles between and the thrust directions. 

5. The tether directions for the minimum thrust-sum configuration are as 

follows. For identical helicopters with equal load sharing, the two tethers must be 
parallel to the sum of the apparent loads owing to the suspended load and spreader 
bar, and very nearly perpendicular to the spreader bar. More generally, if the 
apparent loads are all parallel, the entire system is coplanar with tethers and 
thrust vectors parallel to the common direction; this occurs at hover and is nearly 

obtained at low speeds for any flight condition, helicopter pair, or load distribu- 

tion. At higher speeds, these directions become nonparallel and noncoplanar to 
moderate angles; principally because of the differences between the cargo and heli- 
copter specific aerodynamic forces; however, the tethers are always close to the 
plane of the triangle except near the extreme case of an unloaded tether. 

6. The ratio of thrust requirements for the two helicopters depends solely on 

the spreader-bar tilt, assuming that the tether directions are selected to minimize 
the thrust-sum. A suitable choice for this ratio is given by the helicopters' 
thrust limits, This value is required when the cargo weight is the 

maximum possible for the helicopter pair and is always a possible choice for smaller 
weights. For identical helicopters, equal loading is obtained with the spreader bar 
perpendicular to the apparent suspended load. 

7. If the tether tilt angle ei^ is selected to be different from the value 
for minimum thrust-sum, then the corresponding changes in the equilibrium configura- 
tion can be summarized as follows. First, thrust and tether force vectors vary only 
in their components along the spreader bar; these variations equal in magnitude the 
change in spreader-bar force from its optimum value, but with opposite signs for the 
two tether forces and the two thrust vectors. Second, spreader-bar force and the 
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dependent tilt angles vary linearly with the misalignment (with tan £ 13 ), and 
sensitivity increases with unequal load distribution (with |e£,|)* Third, the 
thrust requirements for each helicopter and their sum, T1, T2, zT, are simulta- 
neously minimized at the optimum configuration in the special case that the 
spreader-bar is perpendicular to ZL, For this case and for small values of 
the sensitivity of thrust requirements to misalignments is much smaller than 
per radian. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California 9^035 
March 3, 1986 
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APPENDIX 


FORCE BALANCE AT THE SPREADER BAR ENDPOINTS 


The forces applied to the spreader bar by the cables (see sketch) are 


F13 = -F13 ul3 
F24 = -F24 u24 
F35 = F35 u35 
F45 = F45 u45 



and the bar's rigid-body equations of motion are 


F13 + F24 + ^ + F45 □ -m^^fb 
0.5 L34 i^ ® (F13 + F35 - F24 - F45) = Hb - Mab 


(A1) 


where Mab is the spreader-bar aerodynamic moment about its midpoint, and ^ is 
its angular momentum rate owing to These equations contain no information 

on the spreader-bar compression, so the force-balance equations for its endpoints 
are written as 


F13 + F35 - F34 i^. = a 
F24 + F45 + F34 i^ = b 


(A2) 


where a and b contain the effects of spreader-bar mass and motion. Appropriate 
expressions for a and b can be obtained by imposing consistency of equation (A2) 
with equation (A1); from which 


a + b = -m^^fb 

0.5 L34 i^ ® (a - b) = Hb - Mab 
and then solving for a, b: 

a = -0.5 m^^fb - i^ ® (^ - Mab)/L34 - 
b = -0.5 m^fb + ij. ® (^ - Mab)/L34 + 


it 


° it 


(A3) 
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where c is arbitrary. However, c can be combined with F34 in equation (A2) 
without loss of generality. 

The results in equation (A3) appear in equations (3.1) in the text in the form 

a = -0.5 = -0.5 mj^(^ + ^) 

b = -0.5 m^jf4 = -0.5 m^j(fb - ^) 

where 

Af = — — i^ ® (Hb - Mab) 

— mj^ -t — 

or, in scalar form, 

'ft = tO’ 

Hb. = E,((|). )[Jb ub. + S(cob. )Jb wb. ] 
t 1 t b b b 

^ ’*’t h * \ 

'^b,N ^ ^2^®t^^3^’*'t^ 

and Afj. can be computed from these equations assuming that Mab^, 4»^, 0(., 9^, 

and are known for equilibrium flight. 

Last, we note that ^ is expected to be negligible (of the order of 0.01 g) 
in equilibrium at all points on a reference trajectory and can be neglected in 
computations. To estimate the size of assume that the spreader bar is a body 

of revolution which is lengthwise symmetric about the midpoint cross-sectional plane 
and has a lengthwise mass distribution that is constant or that declines with dis- 
tance from the midpoint. Then, 


and 


Jb = diag (I^) 

rnaxd.) < m. L34^/12 
i ^ ^ 

IMI ^ ^ (l^l + + 1^) 

\ max / 


(A5) 


where it can be shown for the present context that |ubj “ I'l' |. We expect, for 
reference trajectories within the operational domain and for typical spreader-bar 
lengths and aerodynamics, that equation (A5) yields |^| << 0.1 g. However, ^ is 
retained in the analytical results for completeness. 


75 



REFERENCES 


1. Carter, E. S.: Implication of Heavy Lift Helicopter Size Effect Trends and 

Multilift Options for Filling the Need. Proceedings of the Eighth European 
Rotorcraft Forum, Aix-en-Provence, France, Sept. 1982. 

2. Carter, E. S.; Cooper, D. E.; and Knapp, L. G.: Single Rotor Options for Heavy 

Lift and Potential of Multilift. SAE Paper No. 791087, Dec. 1974. 

3. Sherman, P. F.; et al.: Feasibility Study for Multiple Helicopter Heavy Lift 

Systems. Report R-136, VERTOL Aircraft Corporation, Philadelphia, Pa., Oct. 
1957. 

4. Maciolek, J. R.; et al.: Preliminary Multilift Feasibility Study. Report 

SER 64460, Sikorsky Aircraft Company, Stratford, Conn., May 1969. 

5. Meek, T.; Chesley, G. B.; et al.: Twin Helicopter Lift System Study and Feasi- 

bility Demonstration. Report SER 64323, Sikorsky Aircraft Company, 
Stratford, Conn., Dec. 1970. 

6. Maciolek, J. R.: Master-Slave Controller for Twin Lift Systems. Report 

SER 6647, Sikorsky Aircraft Company, Stratford, Conn., Dec. 1970. 

7. Curtiss, H. C.; and Warburton, F. W. : Stability and Control of the Twin-Lift 

Helicopter System. J. Am. Helicopter Soc., Apr. 1985. (Also, Proceedings 
of the 40th Annual National Forum of the American Helicopter Society, 1984.) 

8. Wilson, G. J.; and Newton, N. N.: Evaluation, Development, and Advantages of 

the Helicopter Tandem Dual Cargo Hook System. AGARDograph CP-121, Feb. 

1973. 

9. Hutto, A. J.: Qualitative Report on Flight Test of a Two-Point External Load 

Suspension System. Paper No. 473» Proceedings of the 26th Annual National 
Forum of the American Helicopter Society, June 1970. 

10. Sheldon, D. F.: An Appreciation of the Dynamic Problems Associated with the 

External Transportation of Loads from a Helicopter — The State of the Art. 
Proceedings of the First European Rotorcraft and Powered Airlift Forum, 
Southhampton, U.K., Sept. 1975. (Also, Vertica, vol. 1, 1977.) 

11. Taylor, M. J. H.: Helicopters of the World. Scribner, New York, 1978. 

12. Davis, J. M.; Bennett, R. L.; and Blankenship, B. L.: Rotorcraft Flight Simu- 

lation with Aeroelastic Rotor and Improved Aerodynamic Representation. 

Vol. 1, TR-74-10A, U.S. Army Air Mobility Research and Development 
Laboratory, Ft. Eustis, Va., June 1974. 


76 



13. Talbot, P. D.; et al.: A Mathematical Model of a Single Main Rotor Helicopter 

for Piloted Simulation. NASA TM-84281, 1982. 

14. Biggers, J. C.; et al.: Wind Tunnel Testa of Two Full Scale Helicopter 

Fuselages. NASA D-1548, 1962. 

15. Shaughnessy, J. D. ; et al.: Development and Validation of a Piloted Simulation 

of a Helicopter and External Sling Load. NASA TP-1285, 1979. 

16. Hilbert, K. B.: A Mathematical Model of the UH-60 Helicopter. NASA TM-8589O, 

1984. 

17. Weber, J. M. ; et al.: A Mathematical Simulation Model of a CH-47B 

Helicopter. NASA TM-84351, 1984. 

18. Houck, J. A.; Gibson, L. H.; and Stelnmetz, G. G.: A Real-Time Digital Compu- 

ter Program for the Simulation of a Single Rotor Helicopter. NASA 
TM X-2872, 1974. 

19. Sheldon, D. F.; and Pryor, J.: A Study on the Stability and Aerodynamic Char- 

acteristics of Particular Military Loads Underslung from a Helicopter. 

Report TN-AM/40, Royal Military College of Science, Shrivenham, U.K., Apr. 

1973. 

20. Pryor, J.; and Sheldon, D. F.: A Wind Tunnel Investigation of Yaw Instabili- 

ties of Box-Shaped Loads Underslung from a Helicopter on a Tandem Suspen- 
sion. Report TN-AM/62, Royal Military College of Science, Shrivenham, U.K., 
Dec. 1974. 

21. Watkins, T. C.; Sinacori, J. B.; and Kesler, D. F. : Stabilization of Exter- 

nally Slung Helicopter Loads. Report TR-74-22, U.S. Army Air Mobility 
Research and Development Laboratory, Fort Eustis, Va. , Aug. 1974. 

22. Matheson, N. : Stability of Helicopter Slung Loads. Aerodynamics Note 364, 

Department of Defense, Aeronautical Research Laboratories, Melbourne, 
Australia, Aug. 1976. 

23. Laub, G. H.; and Kodani, H. M. : Wind Tunnel Investigation of Aerodynamic 

Characteristics of Scale Models of Three Rectangular-Shaped Cargo Con- 
tainers. NASA TM X-62169, 1972. 

24. Laub, G. H.; and Kodani, H. M. : Wind Tunnel Investigation of Aerodynamic 

Characteristics of a Scale Model of a D5 Bulldozer and an Ml 09 Self- 
Propelled 155-mm Howitzer. NASA TM X-62330, 1974. 


77 


1. Report No. 2. Government Accavion No. 

NASA TP- 26 15 

3. Recipient's Catalog No. ' 

1 

4. Title and Subtitle 

General Equilibrium Characteristics of a 
Dual-Lift Helicopter System 

5. Report Date 

July 1986 

6. Performing Organization Code 

7. Author(s) 

L. S. Cicolani and G. Kanning 

8. Performing Organization Report No. 

A-86114 

10. Work Unit No. 

9. Performing Organization Name and Address 

Ames Research Center 
Moffett Field, CA 94035 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Paper 

12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

505-66-01 

[ 15 Supplementary Notes 


Point of contact: L. S, Cicolani, Ames Research Center, M/S 210-3, 

Moffett Field, CA 94035, (415)694-5446 or FTS 464-5446 


16 Abstract 

The equilibrium characteristics of a dual-lift helicopter system are 
examined. The system consists of the cargo attached by cables to the 
endpoints of a spreader bar which is suspended by cables below two heli- 
copters. Results are given for the orientation angles of the suspension 
system and its internal forces, and for the helicopter thrust vector 
requirements under general circumstances, including nonidentical helicopters, 
any accelerating or static equilibrium reference flight condition, any system 
heading relative to the flight direction, and any distribution of the load 
to the two helicopters. Optimum tether angles which minimize the sum of the 
required thrust magnitudes are also determined. 

The analysis does not consider the attitude degrees of freedom of the 
load and helicopters in detail, but assumes that these bodies are stable, 
and that their aerodynamic forces in equilibrium flight can be determined 
independently as functions of the reference trajectory. The ranges of these 
forces for example helicopters and loads are examined and their effects on 
the equilibrium characteristics are given parametrically in the results. 


17. Key Words (Suggested by Author(s)) 

Dual lift equilibrium 
Dual lift 
Slung loads 
Helicopters 


18. Distribution Statement 

Unclassified 

- Unlimited 

Subject Category 01 

19. Security Oasaif. (of this report) 

20. Security Classif . (of this page) 

21. No. of Pages 

22. Price* 

‘ Unclassified 

Unclassified 

88 

. _ 

A03 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1986 



















